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Abstract 



Quantum Chromodynamics (QCD) is the fundamental theory for strong interac- 
tion, one of the four forces in nature. Different from electromagnetic interaction, 
due to its non-Abelian color symmetry, the QCD has the property of asymptotic 
freedom at large momentum transfer, while remains strongly coupled at low ener- 
gies. This property leads to color confinement, i.e. there are no free quarks and 
gluons carrying color degree of freedom, and quarks and gluons are confined inside 
hadrons. In 1974-75, Lee and Collins-Perry suggested that the deconfinement can 
be reached through the ultra-relativistic heavy ion collisions, where the vacuum 
can be excited to a new state of matter or a quark gluon plasma. In the Rela- 
tivistic Heavy Ion Collider (RHIC) at Brookhaven National Lab (BNL), the gold 
nuclei are accelerated and collide head-to-head at the center-of-mass energy of 200 
GeV per nucleon. After the collisions, the huge amount of energy is deposited in 
the central rapidity region to excite vacuum and produce many quarks and gluons. 
In very short time, these quarks and gluons collide to each other and the fireball 
reaches the local thermal equilibrium, and a quark-gluon-plasma (QGP) is then 
formed. After expansion and cooling, the temperature drops substantially and 
the quarks recombine into hadrons which are finally observed in detectors. The 
collective flows such as radial and elliptic flows are observed at RHIC and can be 
well described by ideal fluid dynamics. Further study of the elliptic flow using 
dissipative fluid dynamics gives very low values of the ratio of shear viscosity to 
entropy density, close to the lower bound 1/Att by AdS/CFT correspondence. This 
is one of the most surprising observation at RHIC: the QGP is strongly coupled or 
a nearly perfect liquid instead of a gas-like weakly coupled system by convention. 
The relativistic fluid dynamics provides a useful tool to bridge the gap between 
the initial state of quarks and gluons and the experimental data. Pushed by heavy 
ion collision experiments, the theory of relativistic fluid dynamics also has a lot of 
new development in recent years. This thesis is about the study of three important 
issues in the theory of relativistic fluid dynamics: the stability of dissipative fluid 
dynamics, the AdS/CFT application to shear viscosity in a Bjorken expanding 
fluid, and a consistent description of kinetic equation with triangle anomaly. 

Hydrodynamics is a long-wavelength effective theory for many particle sys- 
tems. The basic hydrodynamic equation consists of conservation equations of 
energy-momentum and charge numbers. One can carry out gradient expansion 
for energy-momentum tensor and charge currents, or equivalently expansion in 
powers of the Knudsen number. The Knudsen number is defined as the ratio of 
the macroscopic length scale (hydrodynamic wave-length) to the microscopic one 
(mean free path). When the macroscopic length scale is much larger than the 
microscopic one, the fluid dynamics is a good effective theory. In expansion, the 
zeroth order corresponds to the ideal fluid. The first order gives the Navier-Stokes 
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equations with dissipative terms like shear and bulk viscous terms. There are a 
lot of candidates for the second order theory. We focus on the widely used Israel- 
Stewart (IS) theory including the simplified and complete version. The simplified 
version can be obtained by correspondence to the macroscopic phenomena. To 
derive the complete version the kinetic theory or the Boltzmann equation is neces- 
sary. In this thesis, we will give the details as how to obtain the complete version. 
The connection between the transport coefficients in the first and the second order 
theories will be demonstrated. 

In the first order theory, causality is violated since it takes no time for a system 
in a non-equilibrium state to reach equilibrium, i.e. the propagating speed of the 
signal is arbitrarily large. In the second order theory, due to finite relaxation 
time introduced, it takes finite time for the system to reach equilibrium. The 
propagating speed of the signal is limited. Therefore the second order theory 
is necessary for causality. However the causality cannot be guaranteed for all 
parameters. The constraints for parameters are then given. We also point out that 
the causality and the stability are inter-correlated. Relativity requires that the 
signal propagate in light-cones. An acausal propagating modes must be forbidden 
by instability or singularity. The connection between causality and stability is also 
discussed. For convenience we work in a general boost frame. It is found that a 
causal system must be stable, but an acausal system in the boost frame at high 
speed must be unstable. 

The transport coefficients can be determined in kinetic theory. There are two 
main techniques to compute the transport coefficients. One is to employ the 
Boltzmann equation, the other is to use the Kubo formula in field theory. We 
will firstly discuss about derivation of the shear viscosity via variational method 
in the Boltzmann equation. Secondly, after a short review of Kubo formula, we 
will compute the shear viscosity via AdS / CFT duality. 

Different from the work of Policastro, Son and Starinets (PSS), we focus on 
strongly coupled QGP (sQGP) with the radial expansion and Bjorken boost invari- 
ance. The information of the sQGP is encoded on the boundary of AdS space via 
the holographic renormalization. Solving the Einstein equation with the boundary 
condition given by the stress tensor yields the metric of the AdS space. In this 
case, the metric is Bjorken boost invariant and has the radial flow. The evolution 
of the shear viscosity as a function of proper time can be obtained via the Kubo 
formula for the retarded Green function given by the AdS/CFT duality. It is found 
that the ratio of the shear viscosity to entropy density is consistent with PSS. 

As another application of AdS/CFT duality, we investigate the property of 
baryons in sQGP in a Wilson-loops-like model. The quarks located at the bound- 
ary of AdS space are connected to a probe D5 brane by superstrings. By studying 
the configurations of baryons with different spins, the screening length of baryons 
can be obtained as a function of spin and temperature. We also study the relation- 
ship between the angular momentum and energy for different kinds of baryons, 
which shows the Regge-like behavior, i.e. the total angular momentum is propor- 
tional to the energy squared. 

As the last topic, we investigate the fluid dynamics with quantum triangle 
anomalies. Generally the relativistic fluid dynamics does not allow the vorticity 
due to parity conservation. Recently it is pointed out that the vorticity has to be 
introduced to relativistic fluid dynamics with anomalies to satisfy the second law 
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of thermodynamics. These new terms are also relevant to the Chiral Magnetic 
Effect (CME) or Chiral Vorticity Effect (CVE). Such terms can be derived from 
the kinetic approach. The coefficients of the vorticity in the case of right-handed 
quarks (or left-handed anti-quarks) and quarks-antiquarks of mixed chirality are 
evaluated. 
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Chapter 1 



Introduction 



1.1 Quantum Chromo dynamics and deconfinement 
phase transition 

1.1.1 Asymptotic freedom 

Quantum Chromodynamics (QCD) is a gauge theory for the strong interaction 
which is one of the four fundamental interactions in the nature. In contrast to 
photons in Quantum Electrodynamics (QED), the interaction for gluons are com- 
plicated because of non-Abelian SU(3) color symmetry. The coupling constant a s 
of renormalized QCD in one-loop approximation is given by 



2'^ 



u s (Q 2 ) = , (l.D 

l + &o^ln(Q 2 /M 2 ) 

where Q is the momentum transfer scale, M is the energy scale and 60 — 33 ^ Nf 
is the first coefficient of /3-function given by the renormalization group equation 
with the quark flavor Nf. Here a s (M 2 ) can be chosen as a s (Mf ) at Z-boson mass 

.1/: mm 

a s {M 2 z ) = 0.1184 ±0.0007 . (1.2) 
The data for a s (Q 2 ) j2] are shown in Fig jl.l| Equation (1.1) indicates that the 



coupling constant decreases with the energy scale Q. This property is the called 
asymptotic freedom |3J H]. The interaction for quarks and gluons will be strong 
in the low energy scale which leads to confinement, i.e. quarks are limited inside 
hadrons in vacuum or the ground state. 

In the early 1970s Lee and Collins et al. [21 [6] proposed that the deconfinement 
can be reached through the ultra-relativistic heavy ion collisions. According to the 
calculations from lattice QCD, the confinement/deconfinement phase transition 
will take place at temperature of about 170 MeV in three flavor case |7] as shown 



in Fig. 1.2 



1.1.2 Experiments for high energy heavy ion collisions 
1.1.2.1 RHIC experiments 

Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) 
has been running since 2000. Two beams of nuclei (Au or Cu) are accelerated to 
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Figure 1.1: Measurements of a s as a function of the momentum transfer Q |2]. 
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Figure 1.2: Phase diagram of QCD from the lattice calculations with non-zero 
chemical potential [8J. 
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Figure 1.3: The elliptic flow measured at RHIC with prediction given by hydro- 
dynamics [llj . 



collide with the center-of-mass energy of 200 GeV/nucleon or 62.4 GeV/nucleon. 
After the most part of nuclei pass through each other, the huge amount of energy 
is deposited in the central rapidity region, and thus excites the quarks and gluons 
from the vacuum. These quarks and gluons form an expanding fireball, and reach 
the local thermal equilibrium within a extremely short time of 1 ~ 2fm/c. It 
is believed that new state of matter, the quark-gluon-plasma (QGP) has been 
formed [HI HO]- The QGP expands and cools down with it freezes out at some 
critical temperature, below which the quarks recombine into hadrons observed by 
the detectors. 

The collective flows such as radial and elliptic flows are observed at RHIC 
and can be well described by ideal fluid dynamics. In non-central collisions, the 
anisotropic momentum leads to the gradient of pressure. This effect can be de- 
tected by analysis of the final particle spectrum in momentum space. The Fourier 
transformation for the particle spectrum in terms of particle azimuthal angle 
with respect to the reaction plane ip r gives 



d 3 N d 2 N iv- l 

E !f = ^^d-y \ 1 + £ 2v " cos[ " W " *' )] (L3) 



where the second coefficient t>2 is the anisotropy parameter, which is also called 
elliptic flow. The data from RHIC have delivered a surprising result that elliptic 
flow v 2 is very large [TT| IT2"t IT3"j Q3] and compatible with the numerical simulations 



of ideal fluid dynamics [T5| ITBl E3 US], see Fig. 1.3 This indicates that the QGP 
is strongly coupled in contrast to the assumption that QGP is a weakly coupled 
system. 

Further study for QGP gives very low values of the ratio of shear viscosity 
to entropy density rj/s ~ 0.1 [1HI [20] close to the theoretical lower bound given 



by AdS/CFT duality |2T], see Fig. 1.4 For a weakly coupled system with well- 
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defined quasi-particles, the ratio rj/s is very large |2"2"| [25] . For the system with 
strong couplings, the lower bound from the AdS / CFT duality is 1/ (47r) [241 125] 
in the comparison with the 1/12 evaluated from the uncertainty principle |26| . 
Therefore, the fact that rj/s of QGP is close to the lower bound is one piece of 
most convincing evidence that QGP is strongly coupled. The calculations for the 
ratio are shown in Chap. |4]for a weakly coupled system and Chap. |5]for a strongly 
coupled system. 

There are a lot of other interesting and important phenomena at RHIC exper- 
iment, e.g. the jet quenching, parton energy loss, heavy quark production, etc. 
(see e.g. Ref. [151 1211 |29j |30j |3TJ [32] for reviews). Those phenomena will also 
be mentioned in the relevant sections. 

1.1.2.2 LHC experiments 

The main goal of the Large Hadron Collider (LHC) experiments at European 
Organization for Nuclear Research (CERN) is to discover the Higgs bosons, su- 
persymmetic particles and other new physics. There are three major experiments: 
A Toroidal LHC Apparatus (ATLAS), Compact Muon Solenoid (CMS), and A 
Large Ion Collider Experiment (ALICE). The purpose of ATLAS and CMS exper- 
iments is to hunt for Higgs boson and new physics, while the ALICE experiment 
is to pin down and study the QGP. 

In the CMS experiment, the two-particle angular correlations for charged parti- 
cles in the proton-proton (pp) collisions at center-of-mass energies of 0.9, 2.36, and 
7 GeV are measured. A long-range, near-side feature in two-particle correlation 
functions have been observed in pp collisions for the first time [33] . A ridge-like 
structure is observed in the two-dimensional correlation function for particle pairs 
with intermediate transverse momentum of 1 — 3GeV/c, 2.0 < |A?y| < 4.8 and 
A(p « with rj the pseudorapidity and <ft the azimuthal angle. Some authors 
[34, 35 J thought that this discovery might imply the QGP has also been formed in 
the pp collisions. The collective flow has also been analyzed [36| |3"7] . The relativis- 
tic fluid dynamics may become a powerful tool to investigate the new phenomena 
in the ultra high energy pp collisions at LHC. 

The first Pb-Pb collision at center-of-mass energy of 2.7 TeV was realized at 
LHC in November 2010 [38l [39] . The ALICE experiment is aimed to search for 
QGP at 2 to 3 times higher temperatures than RHIC. Since the collisional energy 
at LHC is one magnitude larger than at RHIC, the perturbative QCD (pQCD) 
is expected to work much better than at RHIC. The semi-classical Boltzmann 
equation with the collision terms given by the pQCD is a good tool to describe 
non-equilibrium dynamics of the QGO formed in heavy ion collisions at LHC. 

1.2 Relativistic fluid dynamics and kinetic theory 

In long wavelength or small-frequency limit, almost all theories can be described 
by the fluid dynamics as effective theories. L.D. Landau first suggested to apply 
fluid dynamics to the hadronic fireballs [40] . Then Siemens and Rasmussen [H] 
attempted to use the collective transverse flow to describe date of the low energy 
heavy ion collision experiment BEVALAC. Zhirov and Shuryak [42] tried to explain 
the data of the high energy proton-proton (pp) collisions at CERN-ISR using fluid 
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dynamics. Now fluid dynamics becomes a necessary tool to describe data of high 
energy heavy ion collisions [T51 l4*3l HU H5J . 

1.2.1 Second order theory 

The basic equations for fluid dynamics consist of conservation equations of energy- 



momentum and charges ( 2.1|2.2 ). The energy-momentum tensor and the con- 
served currents can be expanded in the terms of the so-called Knudsen number 
pRU W7\ SSI HH], which is defined by the ratio of mean free path to the macroscopic 
characteristic length. The zeroth order of this expansion corresponds to the ideal 



fluid. In the first order, the Navier-Stokes (NS) equations (2.28) are obtained and 
the shear stress tensor and bulk viscous pressure are introduced. The details for 
the expansion in the power series of Knudsen number will be discussed in Sec. 
12331 

There are different representations for the dissipative second order theories of 
the fluid dynamics, e.g. the the theory of the conformal fluid [50], Israel- Stewart 
(IS) theory |4*6| I4"TI 08J H9], the memory function theory [5T| [52] . the extended 
thermodynamics [5TI |53| 154"] . and others [55| 156] . They differ only in non-linear 
second order terms. In this dissortation, we will focus on the Israel-Stewart theory 

only ma unusual]. 

The simplest IS theory is given by a combination of all irreducible quantities 



in the first order theory, see in Sec. 2.2.4 However, this description does not 



demonstrate the fact that quantities in the first and second order theories are 



related to each other. It is not clear whether the simple IS equation (2.32) contains 
all possible quantities in second order theory. The conservation equations can also 
be studied by the kinetic theory (relativistic Boltzmann equation), or the Grad's 

14 moment approximation |l6] (also see |57| 158] in a different metric). A complete 

15 equations are given in a power counting scheme |4Tl US]. The details will be 



shown in Chap. 2.5 The relationship between the transport coefficients in the 



first and second order theory is also discussed in Sec. 2.5.3 



1.2.2 Causality and stability 

It has been pointed out that the first order theory does not obey the causality 



|59| l60| EU |62j [63] . For instance, as will be shown in Eq.(3.4), the heat conduction 
equation in the first order theory is 

at = D w ' (L4) 

where D is the heat conductivity. The dispersion relation in the linear approxi- 
mation is given by 

u = iDk 2 , (1.5) 

which implies that the group velocity of the signal v g = du/dk is proportional to 
the wave-number k. For k — > oo, the group speed goes to infinite and violates 
causality [M]. Therefore, the second order theory is necessary. The second order 



term of d 2 T has to be introduced in the Eq.(1.4) 

o 2 t dT _ d & 2 t 

dt 2 dt dx 



t^ + ^ = D—, (1.6) 
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Figure 1.4: Ratio of r]/s as a function of temperature scaled by transition temper- 
ature T c in different system |20| . 



where T q is the relaxation time. Now the dispersion relation becomes ui ~ ^D/r q k. 
For D < T q , the signal propagating speed is smaller than the speed of light. For 
T q — > 0, the system is acausal again. Therefore, there has to be constraint condition 
for the r q and D, which is called asymptotic causality condition [64, 65, 66j. 

Stability is intimately related to causality [641 165| [66] . A signal propagating 
faster than the light will move out of the light-cone. The acausal propagating 
modes will lead to some non-physical results, e.g. instability and singularities. In 
this case, for all parameters considered the theory will be unstable if it becomes 
acausal. The details will be given in Chap. |3j 

In the linear approximation, the discussion in Chap. [3] will be universal since 
all candidates for second order theories 0H1 SH E01 [SIl E21 E31 [Ml EH El] differ 
only by non-linear second order terms. 



1.2.3 Boltzmann equation 

The transport coefficients can be determined by the microscopic transport theo- 
ries. There are two main techniques to compute these coefficients. The first one is 



to employ the Kubo formula (5.7) where the transport coefficients are expressed by 
the commutator of operators, e.g. the energy-momentum tensor or conserved cur- 
rents [67, 68J . The commutator can be worked out through standard perturbation 
techniques in field theory. 

The second one is to employ the relativistic Boltzmann equation. If the mean 
free path of the particles is much larger than the interaction length, the quasi- 
particle is well-defined. In this case, a semi-classical description for the equation 
of motion of the particle distributions, the Boltzmann equation, works well. As 
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will be shown in Sec. [2X3] and [Bj the energy-momentum tensor and charge 



currents can be determined by the integrals of the distribution functions. By a 
near equilibrium expansion, the transport coefficients can also be related to the 
integrals of the distribution function via Boltzmann equation. 

At high temperature, the shear viscosity in a gauge theory has been found in 
a leading-log form (221 123] 

T 3 

Voc—. — zi , (1-7) 

where g is the coupling constant. For g — >■ 0, 77 — > 00 which indicates that small 
rj/s means a strongly coupled system. Recently, the calculation of shear viscosity 
in a gluon gas with 2 — > 3processes has attracted attention from several authors 
[691 ITOt ITT) . The same technique can also be applied to investigate the transport 
coefficients of dense matter near phase transition [72| 173] . The detail will be 
presented in Chap. |4} 

Along this lines, there are more many other on this topics, e.g. calculation 
of the bulk viscosity \74\ I75| 176] and the second order transport coefficients \77\ . 
There are more applications of kinetic theory to the heavy ion collisions, e.g. the 
collective flow [751 ESI EHl EE], the jet quenching and partons energy loss |H21 
M> ESI ESI EZ1 EFJ EFJ, the shock wave and Mach cone US EH E21 E31 E3 ESI 
and the thermalization (27J EH1 M EH [Ml HD21 [TUB] . 

1.2.4 Numerical simulations of hydrodynamics 

At RHIC experiment, the local thermal equilibrium is established in a very short 
time after collisions. Therefore, the fireball will expand over a sufficient evolution 
time, when the relativistic fluid dynamics works well. In the earlier time of this 
field, the numerical simulations was only for the ideal fluid, see e.g. [T51 03] 
for reviews. However, the ratio of the shear viscosity rj to the entropy density 
s is found to be small but not zero. The simulation for the dissipative fluid is 
necessary. Recently, the simulation for the dissipative fluid dynamics has been 
developed [1041 I105j . In comparison with the t> 2 data at RHIC, the ratio rj/s is 
found to be [106J 

1<4tt-<2.5. (1.8) 

s 

The initial condition is usually given by the Glauber model [15j or Kharzeev-Levin- 
Nardi (KLN) approach [1071 1108"| [T09j. More extended models are also used, e.g. 
MC-Glauber model |110| and fKLN jlllj . Generally, there are mainly three kinds 



of EOS, see Fig. |1.5[ EOS I for the ideal gas of massless partons, EOS H for 
hagedorn resonance gases, EOS Q for a combination of the above two. Recently, 
the SM-EOS Q [HJ UT3], a smooth version of EOS Q, is also used [104J. 

1.3 AdS/CFT 

In quantum field theory, the partition function is written as 

Z QFT = J D(Pe^ d4x( - Co+9Cl) 

= J D<Pe l 5 dixCo [l + ig j (^d + Oig 2 )} , (1.9) 
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Figure 1.5: Equation of state for numerical simulations. 
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where D denotes functional integrals of the fields 0, g is the coupling constant, 
Cq and C\ are Lagrangian for the free and the first order interacting parts, re- 
spectively. For g ^ 1, the theory is weakly coupled, and the perturbation works 
well. If g ^> 1, higher order contributions is not negligible and all orders in the 
expansion should be summed. It is difficult to describe a strongly coupled system 
due to its non-linear and non-perturbative feature. 

In recent years a new technique to deal with strongly coupled systems in gauge 
theory has been developed, made use of the string/gauge duality or the AdS/CFT 
duality proposed by Maldacena and many others |114] I115[ I116| . Here AdS and 
CFT are abbreviations for anti-de Sitter space and conformal field theory, respec- 
tively. In the 't Hooft limit or large N limit (with N the dimension of fundamental 
representation of SU(N) group), the action of the open strings is equivalent to that 
of a conformal SU(N C ) theory. In the classical limit (i.e. the gravitons are almost 
free), the weakly coupled closed strings (free gravitons) in a curved space corre- 
spond to the strongly coupled open strings in a flat space. Finally, the strongly 
coupled field conformal theory in a flat space is equivalent to the theory of closed 
strings in curved space. 

Policastro, Son and Starinets first used AdS / CFT duality to compute the trans- 
port coefficients of sQGP |21j and derived the ratio of rj/s > 1/47T, which is con- 
sistent to the RHIC data. After that, there are more developments in this field. 
Up to date three main applications of AdS / CFT duality have been explored. The 
first one is to use pure AdS/CFT to compute quantities of QCD-like CFT at very 
high temperatures (see e.g. Ref. |117l 11181 1119] for the energy loss, Ref. |120l 1121] 
for the potential of heavy quarks, Ref. |121l 11221 11231 1124] for the screening length 
). The second one is to use the so-called AdS/QCD to calculate the properties 
of hadrons (see e.g. Ref. |125| and Ref. |126| , 1127] for the Sakai-Sugimoto model). 
The third one is the correspondence between gravity and condense matter theory 
(CMT), which is the so-called AdS/CMT, see e.g. Ref. [Ml TM ESS HSU for 
strongly coupled superconductivity and superfluidity, Ref. |132l 1133] for the non- 
Fermi liquids. More reviews on applications of AdS/CFT duality can be found 
in Ref. [131 [1351 EU [137] for the pure AdS/CFT, Ref . [1381 Iffl ED] for super- 
conductivity and superfluidity, Ref. |141] for jets and partons, Ref. [142J for the 
deconfinement phase transition and Ref. [oP] 1143] 1144] for the fluid dynamics. 

In the pure AdS/CFT, two main quantities, the (retarded) Green functions 
and Wilson loops, can be computed via the duality. 



1.3.1 Green functions 

The AdS/CFT duality means that the partition function Zqft[J] i n the conformal 
theory is equivalent to the partition function Z str i ng [(j)\ in the classical string theory. 
In this case, one finds 



JQFT 



[J] = e iSst " n ^ cl] , (1.10) 



where J the source coupled to the operator O and S s t r ing is the action of the 
classical gravity and 



'QFT 



[J] = J D^exp (iS + i j d A x.JO 
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Figure 1.6: Wilson loops in AdS/CFT duality. 



The Green functions in the CFT is associated to the derivatives of the action 
in the AdS space, e.g. the two-point Green function of O is given by the functional 
derivatives of S[4> c {\ with respect to the boundary value of <p, 

5 2 Z QFT [J<P] 



G(x-y) 



t 2 5J(x)5J(y) 

$ 2 S string [0cz] 



J=0 



5J(x)5J(y) 



:i.iD 



z-K), <f>=J 



The problem to compute Green functions in a strongly coupled quantum field 
theory is made to compute the classical action of gravity. 

Based on the work of Ref.[21], we study sQGP with the radial expansion and 
Bjorken boost invariance. The information of the sQGP is encoded on the bound- 
ary of AdS space via the holographic renormalization. The evolution of the shear 
viscosity as a function of proper time can be obtained via the Kubo formula for the 
retarded Green function given by the AdS/CFT duality, see Chap. [5] for detail. 



1.3.2 Wilson loops 

It is well-known that the gauge invariant Wilson loops for quark and anti-quarks 
in QCD can be written in the form 

W[C] = e~ v ^ T = e~ aA ^ , (1.12) 

where C is a contour which is usually chosen as a rectangle in Euclidean space- 
time with the area A(C) = RT, R is the distance between quarks and anti-quarks, 
T is the imaginary time, V is the potential and a = V(R)/R is the string tension. 



As shown in Fig, 1.6 the area of the contour C can be obtained by integrating 
over of the world-volume of the strings. These integrals are given by the Nambu- 
Goto action. In this case, the Wilson loops are found to be 

W[C\ = Z strmg [C] = e"( 5 — , (1.13) 
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where I is the length of the loop C and is the mass field. In order to renormalize 
the potential of heavy quarks, the contributions from the mass of strings l<fi must 
be removed. 

In the standard AdS metric, the potential for quarks and anti-quarks is given 
by [137], 

v q ,{R) - -fffqi — r — > 

where gyM is the coupling constant of the super Yang-Mills theory and also see 
Ref. 1120(1121] for the higher order contributions. This potential is strongly related 
to the phenomena in high energy physic, e.g. the jet quenching |117| 11181 H19j . 

In this thesis, a Wilson-loop-like model will be used to investigate the prop- 
erty of high spin baryons in the QGP in Chap. |5j The quarks located at the 
boundary of AdS space are connected to a probe D5 brane by superstrings. The 
configurations of strings give the screening length of baryons as a function of spin 
and temperature. A Regge-like behavior, i.e. the total angular momentum is 
proportional to the energy squared, is also found. 



1.4 Fluid dynamics with triangle anomalies 

The vorticity vanishes in the first order theory of fluid dynamics due to the parity 



conservation. However, the analysis of the power counting in Eq.(2.83) indicates 
that the vorticity has the same order as other dissipative terms (e.g. shear stress 
tensor or bulk viscous pressure) in the first order theory. It implies that some 



terms related to vorticity are absent in the previous treatment (2.84). 

Recently, the relativistic fluid dynamics corresponding to charged black-branes 
through the AdS/CFT duality was found to have two new terms associated with 
the axial anomalies in the first order theory j!45L 1146] (see, e.g, Ref. |147| about 
the holographic model with multiple/non-Abelian symmetries, or Ref. |148j for 
the Sakai-Sugimoto model). The authors of Ref. |149j have derived these new 
terms in relativistic fluid dynamics with triangle anomalies. Also see Ref. |l50j 
about the similar result obtained in microscopic theory of the superfluid. 

Actually, the anomalous fluid dynamics is closely related to the Chiral Magnetic 
Effect (CME) in heavy ion collisions |151l I152[ I153[ [154J. The separation of the 
left- and right-hand particles or anti-particles leads to macroscopic currents in the 
case of strong electromagnetic field. This effect is called CME. In comparison 
with CME, the separation of left- and right-hand particles given by vorticity is 
also called Chiral Vorticity Effect (CVE) |154 J. The separation of chirality can 
be realized by the change of the topological charge |151| or the quantum triangle 
anomalies |1491 154]. Therefore, it is necessary to introduce an additional term 
proportional to the magnetic field in the conserved currents j M oc B^. And the 
vorticity related to the magnetic field will also be introduced to the conserved 
currents j M oc u^, where 

tut* = e» vaP u v d a up , (1.15) 

is defined in Ref. |149j . Here e fiua ' 3 = 1,-1 for that the order of Lorentz indices 
{lAva0)'\s an even/odd permutation of (0123). 

The anomalous fluid will be studied in a kinetic approach in Chap. |6j The 
transport coefficients for the vorticity will also be obtained. 



Chapter 2 

Basics of hydrodynamics and kinetic 
theory 



Fluid dynamics is an effective theory for any interacting theories in long wavelength 
limit. Recently, it is widely used to describe the data in many aspects of the high 
heavy ion collisions |15 | Ii3 | 1331 135] . 

In this chapter, we will give a brief introduction to the basics of the hydro- 
dynamics and kinetic theory. The basic conservation equations of fluid dynamics 
are given in Sec. 2.1| The first order theory of the fluid, called the relativistic 
Navier-Stokes equations, is shown in Sec. 2.2| Then we introduce the classic and 
quantum Boltzmann equation in Sec. 2J3 and |2.4 respectively. In Sec. 2J3 



we 



give a short review to the complete second order theory of the fluid dynamics via 
kinetic theory. 



2.1 Conservation equations 

The basic fluid dynamic equations are the conservation equations of energy-momentum 
and charge 

= 0, (2.1) 
d»f = 0, (2.2) 

where T^ v is the energy- momentum tensor and j M is the charge current. Here we 
consider the homogenous fluid with only one particle species. 

The tensor decomposition of T^ v with respective to the fluid velocity reads 

= euV - (P + II) A"" + (Wu" + h u u^) + tt"" , (2.3) 

where e, P , II and ix^ u are the energy density, the heat flux current, the 
pressure, the bulk pressure and the shear stress tensor, respectively, and = 
_ U ^ U U is the projector onto the 3-space orthogonal to u^. Here the velocity 
is time-like and normalized to 1, i.e. u^u^ = 1. In a Lorentz boosted frame, the 
velocity is given by 

u" = 7(l,^), (2.4) 
where 7 is the Lorentz factor 7 = yl — v 2 . On the other hand, these thermal 



12 
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quantities can also be expressed by T M ", 



1 

' L 

3 

tt^ = P^Tcf, , W = A" v T va u a , (2.5) 



e = u^T^ , P + n = --AT^ , 



where 

p^afs = 1 ( A m<* a ^ + A ^A m ) - iA M!y A a/3 , (2.6) 

is the symmetric rank-four projection operator. By construction, -k^ is traceless 
7r^ = and ti^ u u u = h^u^ = since m^A^ = 0. 

The tensor decomposition of the conserved current j M reads 

f = nu » + , (2.7) 

where n is the number density and i/* is the diffusion current. By construction, 

n = u,f , u» = A^j v . (2.8) 



Using these tensor decompositions in Eq.(2.1, 2.2), we obtain 

n + nO + fl^i/" = , 

i + (e + P + U)9-d^ -h^-n^d^u, = 0, 
A^ a d a (P + n) - (e + P)u^ - W - A» a h a 

-WB - h a d a u» - A"°0V/9 = , (2.9) 

which represent the conservation of charges and energy, the acceleration of the 
fluid, respectivley. Here A denotes £^A = u^d^A. The 6 = d^u^ is the expansion 
scalar. 

For the sake of simplicity, people usually consider two special frames for the 
fluid dynamics. The first one is the Landau frame or energy frame |155] . In this 
frame, the velocity of the fluid u E is defined to describe the energy flow 

< = , v E = -TW E , (2.10) 
u a E T^T M u E 

which is related to and via 

_L_/i/\ (2.11) 

Then the heat flux current vanishes. The second frame is the Eckart frame 
[156J. In this frame, the velocity is used to describe the charge flow 



with 



u%= ,/" + -!/" 



n 
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Then the diffusion current vanishes. Throughout the thesis, only the Landau 
frame will be used. For convenience, the following quantities is also used 

n 

q» = h» , (2.13) 

e + P v ' 

which is in Eckart frame and — in Landau frame. 

e+P 

For a given fluid (the velocity w M is fixed), there are 15 unknown parameters in 
the equations of fluid dynamics. However, the choice of the frame does not reduce 
the number of the parameters since in Eckart frame the = and in Landau 
frame = 0, then the velocity is not fixed. 



2.2 Navier- Stokes approximation 



2.2.1 Equilibrium state 

The first law of thermodynamics read 



de = Tds + \idn 
dP = sdT + ndfi 



(2.14) 



where s is the entropy density, T is the temperature and /i is the chemical potential. 
The entropy density is given by the Durham-Gibbs relation 



e + P 

s = — — an 

T 



(2.15) 



where a = [i/T. Rewriting Eq.(2.14) with the Durham-Gibbs relation, the follow- 
ing equation is obtained 



d(P(3) = nda — ed(3 
By introducing the new variable /3 M 



(2.16) 



with (3 = T 1 , the Eq. (2.16) becomes 



d(PP") = jgda - T^d/3, 



(2.17) 



where j'q = nu^ and Tq U = eu^u u — P^ u are quantities in the ideal fluid. The 
entropy flow is then 



The differential of the entropy flow reads 

-adj$ + PxdT^ . 



ds% 



(2.18) 



(2.19) 
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2.2.2 Off equilibrium state 



It is straightforward to assume that in an off equilibrium state Eq.(2.19) becomes 

-adf + p x dT x " . (2.20) 



if the system is in a state near the equilibrium one. The entropy flow (2.18) 
becomes 

s » = Pf3» + T Xfl /3 x -af-Q 1 *, (2.21) 

where is the high order deviations — j$ , and T^ v — T^ v . Under the infinites- 
imal changes in a and 0, 

a — » a' = a + 5a , {3\ — > (3' x = (3\ + 5(3\ , 

the change of is 



(2.22) 



Recalling Eqs. (|2.5| 2.8|), we assume that the charge and the energy densities do 

(2.23) 



not change in the off-equilibrium state, 

u,(f-j(;) = u,u u (T^~Tn=0 
but the entropy density flow changes in the way 



. (2.24) 

In principle the temperature and chemical potential will be changed if taking the 
contributions from the second order theory into account, i.e. they are not invariant 
under the transformation of w M . In an off-equilibrium state, the temperature and 
chemical potential are not global, i.e. they are the functions of the location x. 



Taking the additional term into account, Eq. (2.20) and (2.21) become 

s 



n 



and 



dps* 



if - j^a + (T^ A - T^)d^ x - 



(2.25) 
(2.26) 



where g M is given in Eq.(2.13). 



2.2.3 Entropy principle 



Taking Eqs. (2.5, 2.8) in Eq.( |2.26| ) and neglecting Q^, the entropy production rate 
reads 



(2.27) 



The second law of thermodynamics dnS^ > requires z/ M , EE, be in the form 



7T 



fll/ 
NS 



2r]d <>M u u> , 

-CO, 

KTA^d„a 



'ns ~ nj - " ' (2.28) 
which is called the Navier-Stokes (NS) approximation. Here r], (, and k are shear 
viscosity, bulk viscosity and heat conductivity, respectively. Note that in the 
metric convention g^ v = diag{+, —,—,—}, v^d^a must be negative, 

u^a = -nT(Va) 2 . (2.29) 
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2.2.4 Simple Israel- Stewart theory 

In Ref. [46J, the authors suggest that in phenomenology the additional term Q M 
should include all the irreducible quantities in the first order theory 

TQ» = ^((3 U 2 + PiQ 2 + - aoW - a^ x q x , (2.30) 

where and otj are constants with j = 0, 1, 2 and i — 0, 1. 
The entropy production rate becomes [46J, 

V = (3U[-e-P U~2/3 ti-(3oUTd a (/3u a ) + a 1 Td fl (a (3)^ + a d II ^} 
+^[d, L u u - fc-K^ - 2/3 2 7r^ - f3 2 Td a (Pu a )n^ + a^Td^a^u, 
+a x dpv v \ + Pv^-Td^a - - 2/3^ - ^Td^u")^ 

+a 2 UTd fM (a /3) + a ^n + a A T^d v {a x p) + c^tt*] (2.31) 

where the new coefficients dj (i = 1,2,3,4) satisfy a x + a 2 = 1 and a 3 + a 4 = 1. 
Thus the general form for the irreducible quantities are 

C[-0 - 2/3 n + aod^ + aio;ov%] , 
2r\\d <tl u v> - 2(3 2 7r < ^> + a-^u^v^ + ctid^z/^] , 
— /tA^[— T9^a — 2fi]p v + a 2 a ITuj, + a d u Tl 
+a 4 aiiT"u Q + aid a iT"} , (2.32) 

which are the called the simplest Israel-Stewart (IS) equations. Here the approx- 
imation dfj,/3 ~ piifj, has been used, and the coefficients /3j and ccj are assumed 
independent of space-time. The complete IS equations have been discussed in 

Ref. gang @21. 

2.3 Relativistic Boltzmann equation 

The fluid dynamics is closely related to kinetic theory which is based on the Boltz- 
mann equation. The relativistic Boltzmann equation describes the time evolution 
of the single particle distribution function f(x,p), which is based on the following 
assumptions |157j : 

• Only two-particle collisions, or the so-called binary collisions are considered. 

• "Sto/3zahlansatz" collision number ansatz, i.e., number density of binary col- 
lisions at x is proportional to f(x,pi) x /(x,p 2 ). 

• f(x,p) is a smoothly varying function compared to the mean free path l m f p . 

For example, in the 4 interaction the size of the collision region l sc ~ a" 1 / 2 , 
where a is the cross section, needs to be much smaller than the mean free path 
lmfp ~ l/(?2er). In that case, the particles interact in a very small region and 
travel freely in a long distance. For l sc /l m f p <C 1, the semi-classical treatment of 
quasi-particle collisions in Boltzmann equations will be allowed. 
The relativistic Boltzmann equation can be written as 



n = 

^ = 



f = fjV(*,p) = C[f] , 



(2.33) 
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where p^ and E p are the four-momentum and the energy of the particles respec- 
tively, and C[f] is the collision term. With the particle velocity as v p = p/E p , the 
relativistic Boltzmann equation is in the same form as the non-relativistic one 

j t f = d t f + v p -Vf = C[f\. (2.34) 

The collision term C[f] describes the change of the distribution function / in a 
given time dt. For example, considering a two particles scattering 12 — > 3p, the 
collision term is given by 

C[f] = \( dT l2 ^ 3p [f 1 f 2 (l±f 3 )(l±f p )-(l±f 1 )(l±f 2 )f 3 f p ],(2.35) 

* J 122, 

where 

1 3 flU. 

dr 12 ^ P = Wp \T(p, k)| 2 II (2^)3(2^) (2?r)4<y(fcl + h ~ h ~ P) ' (2 ' 36) 
with T(p, k) the amplitude of the scattering. More details on the collision term 



C[f\ will be given in Sec. 4.1 



In the local rest frame, the charge density and the number current are given 

by 

n = j d 3 pf(x,p) , (2.37) 

t = f d 3 p-p*f(x,p) . (2.38) 

In a Lorentz boost frame, the above can be written in a compact form as a Lorentz 
vector 

f = nv?= [ p"f(x,p) , (2.39) 



v 



where j = J , 2 ^~I E , and the four-velocity vector i 



(27r)3£ p 



IS 



1^ = 7(1,^), (2.40) 



with 7 = y/\ — v 2 and v is fluid 3- velocity. 

The energy- momentum tensor can be expressed in terms of /, 

T>» = fp^f(x,p). (2.41) 

The tensor decomposition of T^ v and gives 

n = u^f = J d 3 pf(x,p) , 

e = T^UfM, = J d 3 pE p f(x,p) , 

P = ~^T lw = ~J\t\ 3 f(x f p), (2.42) 
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If the particles are massless, we have the equation of state for the prefect or 
conformal fluid is obtained 

e = 3P . (2.43) 

In relativistic fluid dynamics, the distribution function fo(x,p) in an equilib- 
rium state is 

where a = 0, ±1 are for Boltzmann, Bose and Fermi distributions, respectively 
and 

a - d 9 



(2tt) 3 ' 

with d g the particle's degree of freedom. 
2.3.1 Juttner distribution 

The relativistic Boltzmann distribution is also called the Juttner distribution. For 
the sake of simplicity, all thermal quantities are evaluated in the local rest frame. 
The number density is given by 



n 



A e» /T 4:iim 2 TK 2 (z) , (2.45) 



where K 2 is the Modified Bessel function of the second kind and z = ^. The 
pressure is 

P = A^ /T A-nm 2 T 2 K 2 (z) . (2.46) 
Then the Clapeyron equation PV = nRT is obtained from Eqs.(2.45 2.46[ ) 



P = nT , (2.47) 

where kg = 1. The energy density is 

e = A e fl/T ATT[3m 2 T 2 K 2 (z)+m 3 TK 1 (z)] 
K ( z) 

= 3nT + , (2.48) 
which is identical to the relativistic Boltzmann-Gibbs statistics. 



2.3.2 Conservation laws 

The Boltzmann equation is used to describe a system in or close to an equilibrium 
state. There are conserved quantities because of the time reverse symmetry in the 
microscopic state. Assuming that the quantity \& is conserved during the binary 
collision 12 — > 3p : 

i + ^ 2 = + V P , (2.49) 
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then the macroscopic conservation law will be 



3„ 



\ [ drdPpVtfMi ± / 3 )(i ± /„) - (i ± / x )(i ± f 2 )hf P 

Z J 123 

\ [ dTd 3 P [^l + ^2-^3- *p][/l/ 2 (l ± / 3 )(1 ± fp)} 
4 J 123 



4 

. 



(2.50) 



Choosing ^ = 1, and p 11 , the conservation equations for the charge density and 
energy-momentum are obtained 



d,T^ 



p»d,f(x,p) = 0, 



(2.51) 



The conservation laws follow Eq.(2.49). If Eq.(2.49) is not fulfilled or the time 



reverse symmetry is broken (i.e. the amplitude of the scattering 12 — » 3p is not 
equal to that of the inverse scattering) , the charge and energy-momentum will not 
be conserved. For instance, considering the triangle anomalies (i.e. the right-hand 
quarks will become to left-handed quarks if the topological charge is 1), the time 
reverse symmetry is broken; therefore the quark number is not conserved due to 



drf" = CE B 



(2.52) 



where E and B are the electric and magnetic field, respectively. Here C is the 
constant determined by the quantum field theory. More discussions on this topic 
will be presented in Chap. [6j 



2.4 Kadanoff-Baym equation 

The Boltzmann equation can be derived from the Kadanoff-Baym (KB) equation 
in quantum field theory via gradient expansion. In this section we will derive the 
Kadanoff-Baym equation from Dyson-Schwinger (DS) equation in the closed-time- 
path formalism. Then we will show how the Boltzmann equation can be derived 
from the KB equation. There are many references about the closed-time-path 
formalism and derivation of the Boltzmann equation from the KB equation, see 
e.g. Ref. [T581 



[1591 [1601 rT6Tl fT62l [1681 [16^1 IT651 [1661 fTHTl fT68l H69] . 

We consider a fermionic system in quantum field theory, the two-point Green 
function for a fermion is defined as 



iG{x u x 2 ) = iG{l,2) = (T t [^(l)^ H (2)]) 



(2.53) 



where the operator T t is the causal time-ordering operator and H denotes the 
Heisenberg picture. In order to describe a non-equilibrium system, the closed 
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time path formalism is proposed by Keldysh and Schwinger, see Fig. 2.1 Then 
the two-point Green function on the contour C is written as 

•G(l,2) = (T c [«l)«2)]) 



T c < exp 



i, I d xCi 

c 



^(1)^(2) 



(2.54) 



where / denotes the interaction picture and £j is the interacting part of the 
Lagrangian. The retarded and advanced Green functions are given by 

G R (1, 2) = G(l, 2) - G<(1, 2) , G A (1, 2) ee G(l, 2) - G>(1, 2) , (2.55) 

where Feynman propagator G(l,2) is defined by 

G(l, 2) = fcfa - t 2 )G>(l, 2) + ^ c (t 2 - t 1 )G < (l, 2) 

with 6* c (x) the Heaviside step function in the closed time path C. 



Figure 2.1: Closed time path. 

The Dyson-Schwinger equation reads 

G- 1 (1,2) = G 1 (1,2) + S(1,2), (2.56) 

where Go is the Green function for the free particles and E is the self-energy 
(normally with an i = y/—l as zE) which can be also written as 

E(l, 2) = E 5 (l, 2)5(h - t 2 ) + 0c{ti - t 2 )E > (l, 2) + 6 c (t 2 - ti)E<(l, 2) . (2.57) 

The solution of the Dyson-Schwinger equation can be formally written as 

G(l,2) = G (l,2)- / dA [ d3G (l,4)E(4,3)G(3,2) , 

Jc Jc 

G(l,2) = G (l,2) - f dA [ d3G(l,4)E(4,3)Go(3,2) , (2.58) 

Jc Jc 



Using the Dirac equation and Eq.(2.58), in the case of t\ < t 2) G >,< fulfill the 
following equations 



c/35(l,3)E' 5 (l,3)G < (l,2) 



c 
-ti 



-i 7 ^-m)G < (l,2) 



+ /' 1 rf3[E < (l,3)-E > (l,3)]G < (3,2) 

J t 

+ / >2 rf3E < (l,3)[G > (3,2) -G<(3,2)] , 

J to 

- [ rf3<5(l,3)E <5 (l,3)G < (l,2) 
Jc 

+ /' 1 rf3[G < (l,3)-G > (l,3)]E < (3,2) 

Jt 

+ /" 2 rf3G < (l,3)[E > (3,2) -E<(3,2)] , (2.59) 

J to 
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Note here X 5 is the self-energy depending on a single time corresponding to the 
contribution from the mean field approximation |167] . Without loss of the gener- 
ality, the T, 5 can be neglected in normal situations. 

Taking traces and the Fourier transformation for Eqs.(2.59), we obtain 

Tr[(i 1 ^-m)G < (X,q)] 

^g 1 e^ 1 ^'Tr{[S<(X,g 1 ) - E>(X,q 1 )]G<(X,q)} 



+ 



cV^ie i(<? -^'Tr{£<(X, q)[G>(X, q 2 ) - G<(X, q 2 )}} 

Tr[(-i lfi dg-m)G<{X,q)} 

dy'dq 1 e i ^y'Tv{[G < (X,q 1 ) - G>(X, qi)}E<(X, q)} 
o 



+ 



dy'dqie 



i(q-<n)y' 



Tr{G<(X,q)[Z>(X,q 2 ) -E<(X,q 2 )}} 



where the gradient expansion variables 

Xi + x 2 



X = X 



12 



y = yi2 = x 1 -x 2 , 



have been used. Combining the above two equations, we derive the covariant 
Kadanoff and Baym equation |170| 

z^Tr[7 M G<(X,g)] = Tr[G > S < - G < S > ](X,g) . (2.60) 

It is known that the two-point Green function of free particles must be pro- 
portional to the distribution function f Q of an equilibrium state at the finite tem- 
perature. Therefore, the generalization to the Green function of the interacting 
particles is straightforward. One needs to use the complete distribution function / 
to replace fo in Eq.(2.60). On the left-hand side of KB equation, it will give such 
a term as p^d^f while on the right-hand side, the self-energy S <,;> are associated 
with the amplitude of the interaction. Finally we end up with the Boltzmann 

d 



equation (2.33) 



dt 



f = jfdpf = C[f] . 



2.5 Complete second order theory 

In hydrodynamics, in a small departure from equilibrium, it is assumed that j M 
and T^ y can still provide a complete description of the non-eqilibrium states. Then 
there are 14 independent parameters. Correspondingly, in microscopic distribu- 
tion function, a departure from equilibrium state can also be characterized by 14 
parameters in y, the exponent in the distribution function l/(Exp(y) — a). The 
comparison between the macroscopic and microscopic approach will give part of 
constraints on these parameters. The method is called Grad 14 moments approx- 
imation [46J. 
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2.5.1 First order theory 

In the kinetic theory, the entropy flow is defined as 



where 



/ln^o^-a-^AlnA, 



with 



A = 1 + af, 
y(f) = i;'(f)=ln[A 1 f(x,p)/A(x,p)} 

In an equilibrium state, we have 

Vo = y(fo) = a- p^pP . 



(2.61) 



(2.62) 



(2.63) 



(2.64) 



Since there are 14 independent parameters in j p and T^ u , the derivation from 
equilibrium is given by y — yo which can be decomposed into 



y = (t + a)-{^ + e M )p" + e^p v , 
where e, e M , e pv are 14 small parameters in the first order (with = 0). 



(2.65) 



By using the ra-th and auxiliary moments defined in Eq.(B.l ), the infinitesimal 



changes of j^, T^ u and the 3-rd moment F^ vX with the full distribution f(x,p) 
under an arbitrary variation of f(x,p) are 



5F nuX = F ^fj _ F^ x [f ] = tJ pvX - e p J^ xp + e pa J puXpa . (2.66) 

It is reasonable to assume that the charge and energy density in an off equi- 
librium state should be the same as those in an equilibrium state, 



n 



(u-p)f = n = {u-p)f { 



e = (u-p) 2 f = e = {u-pfh 



Thus, the equation of state will be the same as before 

P(e,n) = P (e ,n ) . 



By construction, Eq.(2.23) becomes 

u^5f = u^uJT^ = 



(2.67) 
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or 



, 2 i a J31J3O — ^41^20 

— \ m + 4— 



D 



A J31J20 ~ J41J10 
4 e* 



20 



D 



(2.68) 



20 



where Eqs.(2.66, B.8) have been used and A* and A** denote A^u^ and A^u^Uy. 

Using Eq.(2.66), all quantities in the first order theory can be written in the 
terms of e M and e^ v . The bulk pressure, shear stress tensor and heat flow are given 
by 



7T 



n 

ft 



5T <nu> = 2J 42€ <xp> 



A p u u^T^ 



e + P 



n 



2 Joi ACry* A^ Q 



where g M is defined in Eq.(2.13) and 

D 31 



A 



j2 ' 



-3 



d In/- 



31 



d In/- 



30 



+ 5 



2^3l(<^21</30 — /20/31) — JAl(J2lJ20 ~ J10J31) g 
/42-D2O 



(2.69) 



(2.70) 



Equivalently, e, e M and e MJy in the Eq.(2.69) can be expressed in terms of II, 7r 
and g M , 



-flU 



fin + e M 
e 



A (3«V - A M ") n - Si« (M g„) + Cqtt^ 

«4 n, 



(2.71) 



where 



A 

A! 

A 2 
Bo 



3A 2 [m 2 + 4 



^31^30 ~~ Jdl J' 



41^20 



D 



20 



-12 
1 



J31J2O ~ <^41 J' 



41^10 



4J 42 fi 
J41 



J: 



#20 

Co 



A 2 



2J42 



31 



AJ ; 



21 



(2.72) 



2.5.2 Second order theory 

For the second order theory (i.e. the theory including the derivative of the quanti- 
ties in the first order theory), one needs to investigate the following 3-rd moment 
with the full distribution function f(x,p) 
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d a F a/lu [f] = I p»p v p a d a f = I p»p v C[f] = . (2.73) 
J p J p 



The contractions of the indices in F atiu and P^ u give 

FIT = [pYf = m 2 f, 
Jp 

P£ = m%f = 0. 
As a rank-2 tensor, P^ v can be decomposed as 

= AuV + BA n" + 2Cu^q v) + Dtt^ 

where A,B,C,D are the integrals of the collision term. 
The rank-2 tensor P^ v can be written in the form 



(2.74) 



pfiu 



(2.75) 



where X^ afS is a rank-4 tensor. Recalling the definition d2.6lh , the entropy prin- 
ciple requires 



< 



ipYd^f = - I e^p"C[f] = e^X^e^ , 



which implies that 



(2.76) 



By the symmetry fl2.76| ) and g^X^ = g aP X^ = 0, the general form of X» ua ? 
is 

X^ = -A ( -Ia^A^ - 3u»u u u a u + A^uV + A a/i vPu v 
3 V 3 



i 



(2.77) 



where A, B, C are the integrals of the collision term C[f] in the Boltzmann equa- 



tion (2.33). 



On the other hand, the tensor d a F a ^ u [f] can be evaluated as 

p^p a f A d a y 
J^ a d a (a + e) - J^d a {(3 p + e p ) + J^ a " a d a e 



pa 



(2.78) 



Inserting Eq.(2.71) and (2.77) into Eq.(2.78) and calculating the derivatives 



of all quantities step by step, the differential equations for the quantities in the 



first order theory give the IS equations which are identical with Eqs. (2.32). The 
details of the calculation can be found in Ref. |4"6| [57] . 
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2.5.3 Complete IS equations 

In the work |16|I57| . the authors neglected the terms a^ v = P^ va ^d a u^^ 9 and d^ot. 
However, it was pointed out in Ref. \£7\ H8] that these terms are in the same order 
as others. 

Generally, there are three length scales in an effective theory for long distance 
limit of a given theory. The first one is the microscopic length scale £ m i cro . In a 
weakly coupled theory with well-defined quasi-particles, this quantity is equal to 
the inter-particle distance. The second one is the mesoscopic length scale £ meS o 
which is identical with the mean free path l m f p in the dilute gas limit. The third 
one is the macroscopic length scale £ m acro ~ e/ \ de\ describing the variety of the 
macroscopic quantities (e.g. the energy density e). Thus, t^ acro is proportional 
to the gradients of the conserved quantities. If the so-called Knudsen number 
K = imeso/^macro is sufficiently small, the expansion in terms of K is equivalent 
to a gradient expansion, i.e., the expansion in terms of powers of l m fpd^- 

The ratios of the quantities in the first order theory to the energy density are 
proportional to K. For example, the ratio of bulk viscous pressure to the energy 
density reads 

U NS (d ■ u C d 

— ~ Imfpdp = K , (2.79) 

where the fundamental relation of thermodynamics, e + P = Ts + fin has been 
used and 

C V 11 11 lmfp , 98n x 

s s \thoT A an\ th X th 
with the help of rj ~ (At^cr) -1 , l m f P ~ (cm) -1 . Here o is the average cross section 



and Xth ~ T" 1 is the thermal wavelength. Note here the result (2.79) is indepen- 
dent on the ratio (/s, i.e., the expansion in terms of K (or the gradient expansion 
) is available in both weakly and strongly coupled theories. The ratios of the cor- 
relations from the second order theory to the energy density are also proportional 
to the K. For example, the ratio of the Tull with the relaxation time Tu for the 
bulk viscous pressure to the energy density reads 

' lliJ n r u u ■ d ~ K^-d ~ K-X th d ~ K 2 , (2.81) 



e e 11 T 4 s 

where the estimation of Tn is employed in Ref. H9] 



0o ~ Yi ■ (2.82) 



It can be proved that the terms <j^ y and 



/IV 



^r^{d a u p -d p u a ) , (2.83) 



give the same contribution as d ■ u ~ <9 M in the power series of K . Therefore, the 
complete IS equations should contain these terms. Taking account of these terms, 
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the complete IS equations is obtained in Ref. |4"?t HH] 

^ = <5-T 7r 7r < ^ > + 2i 7r9 V < ^ > 

+2T 7Tq q < ^u u> + 2T n 7i^u u>x - 2 V 8 2 7r^e 
-2r n n^a u>x - 2K q q < »V u> a + 2A wn n<r^ , 

-Xqqcr^qv + A g nIIV M a + \ qn n^V u a , 
n = U NS - r n n - l Uq d ■ q 

+r Uq q ■ u - c^n^ 

+\ nq q ■ V« + An.vr^V^ , (2.84) 



where 7r^s, IIvs are given by Eq.(2.28) and in the dilute gases limit the 

values of the transport coefficients in the second order theory are given in Ref. |49] . 
These values can also be evaluated by other theories, such as Boltzmann equation 
and AdS/CFT duality. The details for computing the transport coefficients in 
kinetic approach will be shown in Chap. [4] and by the AdS/CFT duality in Chap. 

E 

Note that in dilute gases limit, the ratios of the transport coefficients in the 
second order theory to the quantities in the first order theory are only the function 



of a and f3. For example, as shown in Eq.(2.82), /3q will be determined by the 
macroscopic state , and therefore will be changed with the evolution of the fluid. 
It is a quite different treatment in the simulation for the viscous fluid dynamics, 
for their simplicity, the value of 77 and the relaxation time r n are fixed, see e.g. 
Ref. pun EE]. 



Chapter 3 

Causality and stability 



The first order theory does not obey the causality [5T?| IBTil loTl E21 E3] • The causality 
cannot be satisfied automatically. Therefore, there has to be constraint condition 
for the transport coefficients, which is called asymptotic causality condition |6H 
65, 66j. On the other hand, stability is intimately related to causality [64, 65, 66j. 
A causal theory will be stable. 

In this chapter, we investigate the causality and stability of relativistic dissipa- 
tive fluid dynamics in the absence of conserved charges |65j . In a linear stability 
analysis of the rest frame, we obtain the asymptotic causality condition and find 
that the equations of relativistic dissipative fluid dynamics are always stable. In 
a Lorentz-boosted frame, we find the equations of fluid dynamics are stable if the 
asymptotic causality condition is fulfilled. The group velocity may exceed the 
velocity of light in a certain finite range of wave numbers. However, we demon- 
strate that this does not violate causality, if the asymptotic causality condition 
is fulfilled. Finally, we compute the characteristic velocities and show that they 
remain below the velocity of light if we choose the parameters fulfilled the asymp- 
totic causality condition. The similar discussion for the dissipative currents can 
be found in Ref. |66| . 



3.1 General discussion 

To demonstrate the causal problem in the first order theory, for simplicity only 
the heat conductivity is taken into account in Eq. (B.36) for the ideal fluid. In 
this case, Eq.(B.35) becomes 

= d,f 

= -hoP + hiPO + Jtoa + d^. (3.1) 



Substituting Eq.(3.1) into Eq.(B.33) yields 

1 



a 



D 20 
1 



[(JaiJio-JnJwWe-Jwd^) , 



(3.2) 



D20J20 

where the second line in the massless limit, i.e. J30/J31 = J20/J21 — 3, becomes 

a + -^-[d 2 - (u ■ d) 2 }a = . 
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In the local rest frame, the above equation can be written in the form 

d t a = -T^V 2 a . (3.3) 



If the chemical potential [i is fixed, Eq.(3.3) is actually the heat conduction equa- 
tion 

dT B 2 T . , 

which gives the dispersion relation 

to = iDk 2 . (3.5) 

This implies the system is acausal since that the group speed of the signal is 
proportional to the wave-number k |64j . 

To avoid this problem, the simplest way is to introduce a second order time 



derivative afT in Eq.(3.4). Then Eq.(3.4) becomes a standard wave equation 



d 2 T dT B 2 T 

T «m 2 - + m =D w> (3 ' 6) 



where r q has the dimension of the time. However, there must be some constraint 
condition for r q and D since if r q — > the system will be acausal again. 

3.2 Causality in local rest frame 

The approaches to formulate a second order theory of relativistic fluid dynamics 
is not unique, different approaches differ only by non-linear second order terms 
(see e.g. Ref. HH1 SZlEniEIlESlISnilSllISSlEl])- These differences will vanish 
since a linear analysis is applied here. In that case, the evolution equations of the 
dissipative quantities are given by 

rn-^-n + n = -c<VA 

GST 

t^±» + * = j^V'.. (3.7) 

The investigation of the causality and stability for a hydrostatic background 
in exclusively the low- and high-wave-number limit is given by Hiscock, Lindblom 
and Olson |UU1 [UTj . However, they did not find a generic anomalous behavior of 
the group velocity. The analysis of the causality and stability for the dissipative 
fluid with the bulk viscous pressure only has been done in Ref. |64J, where they 
point out the relation between causality and stability. Since the analysis in the 
fluid with bulk viscosity only is similar to that with shear viscosity only |64j . In 
this chapter, we will make a discussion on the properties of the fluid with shear 
stress tensor only in the absence of the conserved charges. 

For convenience, the following dimensionless parameters will be used 

a = 5, fc.= (e + P)T - = ^, (3.8) 
s n a 
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where e + P = Ts has been used in the absence of conserved charges. A D dimen- 
sional (D > 3) system will be considered. The symmetric rank-four projector is 
in the form 



pHua/3 _ _ ^iia^vP _|_ ^[M0 ^va^ _ A 1 "' A " ' 



(3.9) 



D - 1 

A perturbation around the hydrostatic equilibrium state is introduced 

X = X + 5Xe iult - ikx , (3.10) 

where X = e, 7r Ml/ , with eo = const, tIq U = and = (1, 0, 0, ...). In the linear 
approximation, the perturbation quantities 5X are chosen to be 

S7[ (D-2)(D-2)^ 57r 23 ; S7T 2^ ^ g^D-l) ^ ^ ^^-^(I?-!)^^ 1:L) 

constrained by the normalization condition u^u^ = 1, the traceless condition 7r^ = 
and the orthogonality condition u fJi ir f11 ' = 0. The linearized fluid-dynamical 
equations including the evolution equations (3.7) can be written as 

A8X = , (3.12) 

where the matrix A is in the form 



A 



/TO \ 
5 
G0C0 
\ E J 



(3.13) 



with 





( too 


h o \ 






T = 


-ikc 2 s 

\ o 


h -ik 

r / / 




G 


B = 


diag( J B , . 


• • > -Bo)(D-2)x 


(D-2) , 




C = 


diag(/, . . 


• ) /)(D-3)x(D- 


-3) , 




E = 


diag(/, . . 


• ; /)§(D-2)(D- 


-3)x|(C-2)(D-3) J 






(D-3)x3 



(3.14) 



where c s = \fdP~Jde is the speed of sound and for convenience, the following 
abbreviations have been used 



h = -ik(e + P), 

2(D — 2) 
iu (e + P) , T = -ik — - 1 rj 



f = iu r w + 1 
h 

Ti = —ikr\ 



r, 



D - 1 
2 

V ■ 



D - 1 



(3.15) 
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The determinant of the matrix A must vanish to avoid the trivial solutions of 



Eq.(3.12). Solving det A = gives 



/ 



det T = det 



det B = (det B 
iu /i 
-ike 2 f% — ik 

or/ 



,D-2 



0. 

o, 

. 



(3.16a) 
(3.16b) 

(3.16c) 



These gives the dispersion relation for u(k). 

A nonpropagating mode with the degeneracy (D — 3)[1 + (D — 2)/2] is given 



by the solution of Eq.(3.16a) 



I 

7V 



(3.17) 



The so-called shear modes with the degeneracy 2(D — 2) are given by the solutions 

1 



of Eq.(3.16b) 



2tA V e + P 



(3.18) 



when k is larger than the critical wave number 



le + P Vb 



2r T 



(3.19) 



The numerical result for Eq.(3.18) is shown in Fig. 3.1 



The solutions of Eq.(3.16c) lead to another nonpropagating mode and two 



propagating modes or the sound modes. The analytic solutions in the limit of 
small wavenumber k are 



i 



± k c s + i — k 2 



(3.20) 



while for large wavenumber, we have 



to 



i 



1 + 



± k c s 4 / 1 + 



+ 



1 + 



r c 2 



(3.21) 



where 



2(D - 2) r] 



D-l e+P 



2(D - 2) 
D-l b 



(3.22) 



is the sound attenuation length. The numerical solutions of Eq.(3.16c) are shown 



in Fig. 3.2 The nonpropagating as well as other propagating modes are stable 
around the hydrostatic equilibrium state since all imaginary parts are positive. 
It is identical to the conclusion of Hiscock and Lindblom |59l [60], EI] and others 

|E2E51E5|. 
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Figure 3.1: The real parts (left panel) and the imaginary parts (right panel) 



of solutions for the shear modes obtained from Eq.(3.16b). The parameters are 

47T ' 



a 



6, cl 



| for the 3 



1-dimensional case, D = 4. 
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Figure 3.2: The real parts (left panel) and the imaginary parts (right panel) of the 
solutions for the sound modes (full lines) and the nonpropagating mode (dashed 



line) obtained from Eq.(3.16c). The parameters are a 
3+1-dimensional case, D = 4. 



4vr ' 



6, ci 



for the 
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Figure 3.3: The group velocity for sound modes v g for a = 1/(4-71") , D = 4 , c 2 s = |, 
and 6 = 6 (full line), 6 = 2 (dashed line), as well as 6 = 1.5 (dotted line). If the 



asymptotic causality condition (3.24) is fulfilled, the group velocity will be always 



smaller than the speed of light. Otherwise, the causality is violated. 



3 - 



I - 



k/T 



Figure 3.4: The group velocity of shear modes for D = 4, 6 = 6, c 2 . — |, and 
a = l/(47r) (full line), a = 1/4 (dashed line), as well as a = 1 (dotted line). The 
group velocity is smaller than the speed of light in the large k limit since the 



asymptotic causality condition (3.24) is fullfiled. 



The causality requies that the group velocity 

dReu 



dk 



must be less than the speed of light. For the nonpropagating modes Recu = 0, 
the causality is associated with the behavior of the imaginary part [M]- Generally 
speaking, a k 2 dependence of any nonpropagating mode can be considered to 
violate the causality. In that case, the results in Eq. fl3~T7|3.20[3.2ip show that 
two nonpropagating modes are causal. For the sound modes, the group velocity is 



shown in Fig. 3.3 The group velocity for the shear modes (3.18) is shown in Fig. 
3.4 Unfortunately, there are divergences near the critical wave number k c . The 



details for these divergences will be discussed in Sec. 3.4 



Taking solutions for all modes into account, it is found that the causal condition 
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is determined by the behavior of group velocity in large k limit, 

< = <^ = >™ ^ = C . (3.23) 



Consequently, the asymptotic causality condition reads 

For conformal fluids, the above causality condition is always satisfied since 

b = 2(2 - In 2) ~ 2.614 > 2 , (3.25) 

where c 2 s = 1/(1?- 1) and the values of a = rj/s ~ l/(47r) and t w = (2 -In 2)/(2irT) 
are derived from the AdS/CFT correspondence [50| I171[ I172j . 

The stability of a propagating mode is associated with the the behavior of 
imaginary part. If the imaginary part of the frequency in this mode is always 
positive, the damping amplitude of perturbations will decrease with the evolution 
and the system is therefore stable. 

Note that in local rest frame the stability of relativistic dissipative fluid dy- 
namics is not affected by the causality. The stable propagating modes with acausal 



parameters can be observed in Figs. 3.1, 3.2 For instance, an acausal fluid with 



D = 4 and b = 1 is demonstrated to be stable for the shear modes in Fig. 3^ and 
sound modes in Fig. 3.2[ respectively. 



3.3 Stability in Lorentz-boost frame 

In order to demonstrate the intimate relation between causality and stability, it is 
necessary to consider the system in a moving frame. For simplicity, the space-time 
dimension is restricted to be D = 4. 



3.3.1 Boost along x direction 

The perturbation of the fluid velocity in a frame boosted with the constant speed 
V along the x direction is given by 



u 



I fx iujt—ikx 



(3.26) 



where 



u 



1 A' 



7v (l, V, 0, 0) , 5v! " = (VyvSu*, iybu x , 5u y , 5u z ) . (3.27) 



The linearized fluid-dynamical equations are given by AX = with 

X = (5e, 5u x , 5tt xx , 5u y , 5-K xy , 5u z , 5n xz , 5ii yy , 5-K yz ) T , 

and 



.4 



( T x \ 

B x 

Gi C x 

\ Ex J 



(3.28) 



(3.29) 
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where 
T x 
B L 

C x 
with 




T 12 



prffv{^V — k) 



diag(S i,-B i 
2 





Ei = F 



irjjviuV — k) 



i-f v 2 V(coV - k) 
ijy 2 (uV — k) 

lv 2 F 
i lv {uj-kV){e + P) 
ir\^y{ujV — k) 







i(uV - k) 
F 



(3.30) 



T n 

T 12 

Ti 3 

T u 
F 



+ V z c z s ) - ikV(l + c 2 SJ . 



-2.2\ A foir, i - 

i[2uV -k(l + V 2 )](e + P) , 
iuV(l + c 2 s )-ik(V 2 + c 2 s ) , 
i [w(l + I/ 2 ) - 2A;7] (e + P) 
i7y(o; — £7)7"^ + 1 . 



(3.31) 



The determinant of the matrix A is detA = det7\ x detl?i x F 2 . From F — 0, 
two trivial propagating modes are found 



to 



+ kV 



(3.32) 



which correspond to nonpropagating modes in the local rest frame. From det B\ 
0, four modes corresponding to the shear modes are observed 



iT - 2o(l - b)W lv ± y'-T 2 + AiakTV lv l + Aa 2 bk 2 -f v 2 

(3.33) 



2a(b - V 2 ) lv 

On the other hand, the sound modes are given by the following equation 



c 2 s (e + P) [1 - HvrAkV - u)] {k 2 [V 2 + (7 - 1) V7 2 , + l] 
+ 2kVu [(V - 1)7 7 2 - 1] + VW - c- 2 {u - W) 2 } 



+ -i lv r){k - Vuf {kV [c 2 sl 2 v V(l -V)-l]+u} = 0. 



(3.34) 



In contrast to the results in the local rest frame, the appearance of negative 
imaginary parts (i.e. the theory becomes unstable) could be observed in the right 



panel of Fig. 3.5 where the parameter set violates the asymptotic causality con- 



dition (3.24) 



3.3.2 Boost along y direction 

The velocity in a frame boosted with the speed 7 along y direction is given by 



in 



7v(l,0, 7,0) , 5u'v = iy lv 5u y ,5u x , lv 5u\5u z 



(3.35) 
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Figure 3.5: The imaginary parts of oj for a boost in x direction with velocity 
V = 0.9 . In the upper panel, the parameter set a = b = 6, (? s = | fulfills the 



asymptotic causality condition (3.24), while in the lower panel a = 
| violates this condition. The dashed lines are for the shear modes 
lines are for the sound modes. 
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while the solid 



CHAPTER 3. CAUSALITY AND STABILITY 



37 



The matrix A is 



with 



and 



/ T 2 

G 2 



#6 

H 7 



H 2 

C 2 




\ 


E 2 J 



(3.36) 



T 9 



H, = 



1^(1 + c 2 s V 2 ) 
— i/cCg 


2iu;y( £ + P) 7 ^ 


— \iuVri^ v 



-ik lv (e + P) 
ioj^yie + P) —ik 
— ^ikr) Fi 

-ikV 
iuV 




B 2 
H 7 



( luV^il + c 2 ) -ikV-yv(e + P) \ 

^7^77 



V / 

/ ico^(l + V 2 )(e + P) -ik 

—ik'jvV Fi 

iu-i V {e + P) 

y -ikr] 

( |iwy 7 ^ ooo), 

( iuVrtv ) , 



o \ 


—ik 
F 1 ) 



(3.37) 



H 2 = ( iuV 2 ) T , 

H 4 =(iuV ) T , 

H 5 = (0 ) T , 

G 2 = ( 1^7^ ) , 
C 2 = E 2 = F 1 , 

where F\ = iwyv^n + 1- 

Solving the equation det A = gives all modes. The nonpropagating 
takes the same form as in the local rest frame, 



u = 



i T ± V-T 2 + 4a 2 bk 2 - Aa 2 k 2 V 2 



The shear modes are given by the solutions 

1 

U± ~ 2a(b - V 2 ) lv 
and the equations for the bulk modes are 
3c 2 s (e + P)(-i + lv r^){k 2 + VW v uj 2 ) 

+ 7yw {Ak 2 r] + 7y w [3i(e + P) + AV 2 -f V r]UJ - 3(e + P)7 y r^] } = , 



(3.38) 
mode 

(3.39) 
(3.40) 



(3.41) 
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respectively. 

The numerical results are shown in Fig. 3^ with the conclusion same as the 
case of boost along x direction. It is confirmed that if the asymptotic causality 
condition (3.24) is fulfilled, the system will be causal and stable. However, the 
reverse is not true. A stable theory may also violate the asymptotic causality 



condition (3.24) 




k/T 



Figure 3.6: The real and imaginary parts for the shear modes (dashed lines) and 
sound modes (solid lines), for a Lorentz boost in y direction with a = b 



6, c 2 s = l,V = 0.9,D = 4. 



3.4 Divergences in shear modes 



As mentioned in Sec. |3.2| (see also Fig. 3.4), there are divergences near the critical 



wave number k c . However, the analysis for the stability of the fluid in a moving 
frame implies that these divergences do not affect the stability. Moreover, a lot 
of work [641 165| 166] show that the causality of theory is guaranteed if the group 
velocity in large k limit is smaller than the speed of light. This problem has also 
been studied in the classical electrodynamics (e.g. see |173j ). The divergent group 
velocity may become superluminal when one analyzes the propagating modes of 
electromagnetic waves in some special material. However, such kind of divergences 
is considered to be unphysical. 

The similar analysis for the divergent group velocity indicates that the diver- 
gences in shear modes do not affect the causality of the theory. The perturbation 
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SX in Eq.(3.10) is given by 



duo SXj(u) 



(3.42) 



where j denotes the different modes and kj(u) is the inverted u(k) of the respective 
mode. The inverse Fourier transform gives the components 5Xj(uo) 



,uo 



1 

2^ 



dtSX{0,t) e 



-iivt 



(3.43) 



where 5X(0, t) = for t < as a result of the assumption that there is no change in 
a fluid-dynamical variable before t = 0. It is found that the • SXj(uo) is analytic 
in the lower half of the complex uo plane. If the asymptotic causality condition 
is satisfied, the imaginary part of uo is positive; therefore the singularities only 
arise in the upper half-plane and the system will be stable. On the other hand, if 
the asymptotic causality condition is violated, the singularities may appear in the 
lower half-plane and the system is unstable. 

In order to demonstrate that the divergences of group velocity in shear modes 



do not violate the causality, it is necessary to compute the contour integrals (3.42) 
in the complex uo plane. To close the contour, the asymptotic behavior of the 
dispersion relations, i.e., the behavior in large k limit, must be known. In large k 



limit, the exponential in Eq.(3.42) becomes 



with 



exp[iuot — ikj{uo)x\ — > exp 



lim KeUj(k) 

k— >oo 



UO 



ill 



[X 



V™ k 



(3.44) 



(3.45) 



where is given by Eq. (3.21) 



a.) 



If x > Vg S jt, the integral contour must be closed in the lower half-plane. If 



the asymptotic causality condition is fulfilled, Eq.(3.42) vanishes since there are 



no singularities in the lower half-plane. If x < t, the integral contour must be 



closed in the upper half-plane. Therefore the value of SX(x,t) in Eq.(3.42) will 
be nonzero because of the singularities. On the condition that the asymptotic 
causality condition is fulfilled , i.e., the asymptotic group velocity is smaller 
than the speed of light, the locations x lie within the light cone. In that case, the 
system is causal. 

The conclusion is that the causality of theory as a whole is guaranteed by the 



asymptotic causality condition (3.24) 



3.5 Characteristic velocity 

The fluid-dynamical equations with nonlinear effect can be written as 

(Ald t + Al b d x + Al b d y ) Y b = B a , (3.46) 
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with y b T 



TT xy ) and I?J = (0, 0, 0, tt xx , ir xy ). The expressions of 
the matrix A are given in the Appendix |XJ The characteristic velocities are given 
by the solution of the following equations |59l 1601 1611 1621 163] 



det(v x A*-A x ) = 0, 
det(« y A* - A") = 0. 

In the local rest frame, the characteristic velocities are given by 

' , 



(3.47) 



± 
± 



(3.48) 



The numerical results for the b dependence of one of the five characteristic veloc- 
ities are shown in Fig. |3.7 



The characteristic velocity which includes all the non-linear effect of the fluid 
dynamics also show that the system will be causal if the asymptotic causality 
condition is fulfilled. 



1.5 




0.5 - 




Figure 3.7: One of the five characteristic velocities given by Eq.(3.47) with w M 
(V5/2, 1/2, 0), tt xx = n xy = 0, and c 2 s = 1/2. 
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3.6 Discussion 

In this chapter, the analysis of causality and stability of the fluid dynamical equa- 
tions is performed. Considering a linear analysis, the so-called asymptotic causal- 



ity condition (3.24) is obtained. The divergences in shear modes are also observed. 
However, the analysis from the contour integrals show that the causality as a whole 
is determined by the asymptotic causality condition. The stability is found to be 
intimately related with the causality of the system in a Lorentz boosted frame. 



The system will be always stable if the asymptotic causality condition (3.24) is 
fulfilled. More work on this topic could be found in Ref. [65J for the competition 
of bulk and shear viscosity, Ref. |B6j for the analysis of the fluid dynamics with 
the heat conductivity only. 

From the asymptotic causality condition, it is found that the NS equation is 
acausal if the relaxation time for shear viscosity goes to zero. On the other hand, 
it also implies that the equations of second order theory are not automatically 
causal by construction. It is easy to check that the results from the Grad's 14 
moments HH] as well as the results from the AdS / CFT fulfill the asymptotic 



causality condition (3.24) 



Chapter 4 

Transport coefficients by Boltzmann 
equation 



In chapter 2, we have introduced IS theory of hydrodynamics [4"6l l4"9l [57] . There 
are many transport coefficients in the theory. These coefficients cannot be deter- 
mined in hydrodynamics but can only be determined in the underlying microscopic 
theory. The collision term in the Boltzmann equation has the microscopic nature 
since it is given by the invariant amplitudes of microscopic processes. In this 
chapter we will discuss about the procedure of computing these coefficients in the 
Boltzmann approach based on Ref. [7Q"| I7T| I73| 172] . 

In a leading order expansion of the coupling constant, there are an infinite 
number of diagrams in non-equilibrium field theory |174[ 1175] . However, it is 
proven that the summation of the leading order diagrams in a weakly coupled 
4 theory pH HZQ Q2B1 EH EE] or in hot QED [179] is equivalent to solving 
the linearized Boltzmann equation with temperature-dependent particle masses 
and scattering amplitudes. The conclusion is expected to hold in weakly coupled 
systems and can as well be used to compute the leading order transport coefficients 
in QCD-like theories [221 EDI GU [M USD], hadronic gases [Ml [Ml [Ml [Ml [Ml 
□EE]- 



4.1 Order expansion 



The basic feature of relativistic Boltzmann equation (2.33) is shown in Sec. 2.3 
In some case, we can add the external field F^ v = d^A u — d u A 11 to the Boltzmann 
equation 

which can be written in a classical form. Here q is the conserved charge, qF^^- = 
g(E+v' x B) is the classical Lorentz force with v' = p/E p the velocity of a particle. 



The distribution function fo(x,p) in an equilibrium state is given in Eq.(2.44). 
It is straightforward to extend it to an off-equilibrium state is 

f( x >P) = p/ 3(u-p- M )+ x (x,p) _ ~ ' ( 4 - 2 ) 
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where x ~ ®n u v ~ 0(K) is a small quantity. The complete distribution function 
f(x,p) can be expanded in terms of the power series of K around an equilibrium 

state 



f(x,p) = f (x,p) + f 1 (x,p) + O(K 2 ) 



(4.3) 



where 



h = -X(x,p)f (x,p)[l + af (x,p)] , (4.4) 
Substituting Eq.(4.3) into the Boltzmann equation ( 4.1[ ) yields the Boltzmann 

dfi 



equation in the i-th order 



dt 



am 



where 



Q[f] = c[n) !/*-«] . 



(4.5) 



(4.6) 



n=l 



with C^ n ' the n-th order derivative of the collision term C[f]. Note that the energy 
conservation in the process of two particles scattering ab — >■ cd, i.e. E a + E b = 
E c + Ed leads to the vanishing of C[/o], 



C[f ] = \f rfr a ^ cp [/ V fe (l + a/o c )(l + a/o P )-(l + a/o)(l + a/o)/o7 ( 
= ^ / dT ab ^ cp (l + afc)(l + afP)(l + af«)(l + af b ) 

1 Jabc 



X 



fa fb 
JO JO 



f C f P 



;i + af$)(l + af b ) (l + a/ c )(l + a/ p ) 



. 



Here we have used 



fa fb 

Jo Jo 



fc fP 

JoJo 



l + afS)(l + af b ) (1 + afS)(l + afg) 



(4.7) 



The derivative of the fo(x,p) reads 



dfo 
dt 



-f«[l + ar*(x,p)) \ ^d, [/3(u ■ p - nq)] 



+ [g(E + VxB)]-V p {P{u-p- m )\ 



It is convenient to evaluate the above equation in the local rest frame (i.e. u° — 1, 
u = 0), 



dfo 
dt 



-fe«[l + af e «(x,p)} 



LRF 



p % pp 



E p d t f5 + p l (di/3 + f3d tUl ) + ^-$d L u 



E,. 



p^ 
E~, 



-qd^a + /%(E + v' x B)] 



E„ 



(4.9) 



Using the Eq.(B.41) in the local rest frame, and the following identity 
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P yd jUi = Ly - ^p 2 J a, + ^P 2 V • u 

with (Tij = \ (djUi + diUj) — 4<%V ■ u, we obtain 



(4.10) 



dfo 
dt 



-/o[l + af (x,p)} 



LRF 



n q 
e + P + K 



p • Va + p 



E 

T 



(4.11) 



where the particles are assumed to be massless and therefore the bulk viscous 
pressure is vanished. 



Because of the Boltzmann equation (4.5) the dissipative term \ hi Eq.(4.2) can 
be written as 



X 



1 

T 



n q 
+ 



E 



r + P ' EJ V P" Va + P T 



(4.12) 



which indicates the fact that the heat conductivity is equal to the electric conduc- 



tivity. By considering the contribution from the external fields, the in Eq.(2.28) 
now becomes 

= l ;Ty i " [d v a+^j . (4.13) 



4.2 Shear viscosity 

To show the details of computing transport coefficients via the Boltzmann equa- 
tion, we choose the shear viscosity as an example. For a review of the shear 
viscosity, see e.g. Ref. [187J. 

4.2.1 Non-relativistic system 

Suppose the fluid is following along the y direction with fluid velocity v y (x) which 
is a function of the transverse position x. The friction force per unit area felt in 
the yz plane is proportional to the gradient of v y along x, 



A 



r]d x v : 



v ■ 



(4.14) 



where rj is the shear viscosity and its inverse is called the fluidity. In the molecular 
theory of dilute gases, one can estimate the value of shear viscous coefficient rj as 
Maxwell did. The number of particles which are moving through the unit area in 
yz plane in unit time is 

AiV = v x fdv x dv y dv z , (4.15) 

with / the distribution function. The total momentum transferred in unit time 
gives the force 

F r°° 

(4.16) 



F 
A 



mVyAN 
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which vanishes if using the equilibrium Maxwell distribution function. It implies 
that the shear viscous effect is a dissipative phenomenon in an off-equilibrium state. 
In order to evaluate this effect, an off-equilibrium distribution function / = /0 + /1, 
where fo is the equilibrium one and f\ is the fluctuation near the equilibrium state, 
is considered. The linearized Boltzmann equation d t f + v • df = gives 

v x d x f = --h , (4.17) 

where the formula in the right hand side is given by the assumption dt(f — f°) = 
— (/ — /°)/to with To the time in which the system comes back to the equilibrium 



state. The solution of Eq.(4.17) is 



df dvy , . 

/1 = -T v x d x f = -r v x . (4.18 

oVy ax 



Substituting the above solution into the expression (4.14) yields the shear viscosity 



77 = — j mvlv y T *^-dv x dv y dv z = mr J vlf dv . (4.19) 

More simplification will be taken. In non-relativistic statistic physics, it is known 
that 

m j vlf dv = nmvl = nk B T , 

where the k B is the Boltzmann constant and n is the particle number density. 
Thus, the shear viscosity will be 

v = r nk B T ~ l -^nk B T oc Vf , (4.20) 

v 

where l m f p is the mean free path which is proportion to n" 1 . This result indicates 
that in non-relativistic case the shear viscosity is only a function of temperature 
T and is independent on the number density n. By using the approximation that 
v x ~ |u and v x = v x , the famous formula given by Maxwell is obtained 

V = T^nplmfp ■ (4.21) 

Instead of the picture of quasi-particles, Frenkel and etc. |188| gave a simple 
picture for the motion of liquid molecules. The shear viscosity in a liquid is given 
by 

77 ~ hne E,{kBT) , (4.22) 

where h is the Planck constant and the collision time of the molecules was as- 
sumed to be h/{ksT) which is the shortest timescale in the liquid. In contrast 



to the results (4.21), the shear viscosity in the liquid decreases with temperature. 
Therefore the value of the shear viscosity must be minimum in the critical point 
of liquid-gas phase transition. 



From Eq.(4.21) and Eq.(4.22), the ratio rj/n (or 77/p, the kinetic viscosity with 
p = run) is found to be a good quantity to describe the minimum value in the 
critical point. In QGP created by heavy ion collisions only net number of quarks 
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is well-defined. The exact number of quarks and gluons is unknown. In that case, 
the ratio r)/n is not good enough to describe the properties of the fireball. The 
entropy density s is well-defined and proportional to the n. Therefore, one can 
choose rj/s around the phase transition to replace rj/n. The uncertainty relation 
gives pl m fp — h, thus, rj/ s > h/ks ~ 1 [26] which is an estimate of l/(47r) from 
the AdS/CFT correspondence [2T1I25]. 

4.2.2 Relativistic system and variational approach 

In this section we will introduce the variational method in the Boltzmann approach 
to shear viscosity. We now turn to a relativistic system. A good example is the 
shear viscosity for a quark gluon system [701 I7T| [721 [75] . 



The shear viscous term in Eq.(4.12) can be rewritten as 



where the heat and electric conductivities are ignored and 



with 



Bij = A x {p) \pf Iij , 



(4.23) 



(4.24) 



(4.25) 



series of K yields 



Substituting Eq.(|4.23) into the expression of T^ u up to the first order of the power 
V= l^fl E~J Ap |P|2 IijBij{p) = B) ' (426) 



where S and B are matrices with the components Su = |p| Jj,-. 



On the other hand, the linearized Boltzmann equation (4.5) with x i n Eq.(4.23) 
reads 



I I 2 T 

P h 



^ f rfr 12 , 3p A 1 A 2 / 3 (Ap)- 1 [B^p) + B tJ (k 3 ) - B l3 (k 2 ) - Bijih)} , 

^ J 123 

(4.27) 

which involves collision terms. The above equation can be written in a compact 
form 

S = CB , (4.28) 



where the matrix C is determined by Eq.(4.27) 



By using Eq.(4.28), Eq.(4.26) becomes 



20T 



dr 12i3p AiA 2 / 3 / P [Bijip) + B i3 (k 3 ) - Bij(k 2 ) - B^h)} B^p) 



123p 



' ' dT Ut3p A 1 A 2 f 3 f p [Bijip) + Btjiks) - Byfa) - B^h)} 2 



8°^ Jl23p 

(B,CB) 



(4.29) 
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which implies that the matrix C is positive. Then we obtain 

(S, B) = {B, CB) . (4.30) 

A straightforward way of computing the shear viscosity is to employ the so- 
lution to Eq. (4.28) in Eq.(4.26). However, there are two problems. The first 
one is from the numerical technique. The numerical errors in this kinds of inte- 
gration equations will lead to some kinds of the divergent behavior |2"2"| |2"3"] . The 



second problem is the solutions of Eq.(4.28) might not fulfill the Eq.(4.30) since 
the solutions of integration equations are not unique. 

Instead of solving Eq.(4.30) directly, the variational method is normally used. 
Eq.(4.29) can be rewritten as 

V 



2(S,B)-(B,CB) 
(S, C~ l S) - (A,CA) 



(4.31) 



where A = B - C^S. If Eq. 04.27) is not fulfilled, r] < (S^^S) because of 
positive C. In numerical calculations, one needs to obtain the maximum value of 
Eq.fl4~31~| 1221123]. 



In variational approach, we solve Eq. (4.30) instead of Eq. (4.28). The critical 
step is to find a good form of Ai(p) to make 77 as large as possible (see e.g. 
[701 EJ H2J US I182j ). As an assumption, one could expand Ai{p) by a set of 
orthogonal polynomials 



A 1 (p) = \p\ y ^Tb r B^\f3p). 



(4.32) 



r=0 



where B^ r \(3p) is a polynomial up to {(3p) T and b r is its coefficient. Here y is a 
constant to make the numerical error get the fastest convergence [72 ^ I182| . In the 
numerical calculations, y is chosen to be 1 or 2. The orthogonal condition is set 
to be 



15T 



IP 



2+y 



f p A p B( r \pp)B( s \pp) = B^8 r , s 



(4.33) 



where B^ is constant depending on the integrals in Eq. (4.33). W ithou t loss of 
generality, we can assume BM = 1. Substituting Eq.fl4.32k into Eq.(|4.29|) yields 



Vtest = ^ b r b s C rs =< b\C\b > 



(4.34) 



where \b >= (bo,bi, b rmax ) T , the inner product is defined as < A\B >= AB and 
1 



80T 

4 



rfri2,3pAiA 2 /3/p 



123p 



X 



(-l) m+n IiAPm-Pn)(\Pn\\p m \) y B^\^p m )B^(Pp m ) 



m,n=l 



where Pi = (p,ki,k 2 ,kA and C is found to be a positive constant matrix into 

(4.35) 



Eq.(4.26). Inserting Eq.(4.32), we have 

J2 b rS r =< S\b > 



Vtest 
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where \S >= (So, Si, S rmax ) T with the r-th component 

S r=Y^ff W fpAp \P\ 2+VB(r) ^P) ■ ( 4 - 36 ) 



From Eq.(4.34) and (4.35), \b >= C'^S >, we obtain 



Vtest = <S\CT 1 \S> 

2 



E 



„ K 15TE P 



' 1 f p A p ) \p\ 4+2y BV(Pp)BW(Pp)C- l (p) 



According to the orthogonality condition (4.33), onlyl^ ) will survive in Eq.(4.35) 
then we finally obtain 

Vtest = J ^fY p fAp IP|2+% = b ° Bi0) = W^ 5 " 1 )" • ( 437 ) 

Moreover, it is proved by the authors of Ref. [721 173] that the value of rjtest 
will increase with r max increasing. Therefore, the value of r max depends on the 
numerical precision. 

4.3 An example: shear viscosity of a gluon plasma 

Recently perturbative QCD calculation of rj/s of a gluon plasma has raised wide 
attention. Xu and Greiner (XG) used a parton cascade model to calculate rj/s 
[69| 1189] . They claimed that the dominant contribution comes from the inelastic 
99 ^ 999 (23) process instead of the elastic gg — > gg (22) process: the 23 process 
is 7 times more important than 22. This result is in sharp contrast to AMY's 
result [221 123] where the 23 process only gives ~ 10% correction to the 22 process. 

Both XG and AMY use kinetic theory for their calculations. The main differ- 
ences are [70] [71] (i) XG uses a parton cascade model [190 J to solve the Boltzmann 
equation and, for technical reasons, gluons are treated as a classical gas instead 
of a bosonic gas. On the other hand, AMY solves the Boltzmann equation for a 
bosonic gas. (ii) AMY approximates the Ng <-> (N+l)g processes, N = 2, 3, 4 . . ., 
by the g -H- gg splitting in the collinear limit where the two gluon splitting angle is 
higher order. XG uses the soft gluon bremsstrahlung limit where one of the gluon 
momenta in the final state of gg — > ggg is soft but it can have a large splitting 
angle with its mother gluon. 

In an earlier attempt to resolve the discrepancy between XG's and AMY's 
results [70], a Boltzmann equation computation of 77 is carried out without tak- 
ing the classical gluon approximation (like AMY's approach) but the soft gluon 
bremsstrahlung limit is applied to the 23 matrix element. It was found that the 
classical gas approximation does not cause a significant error in 77/s. However, the 
result is sensitive to whether the soft gluon bremsstrahlung limit is imposed on the 
phase space or not. If this limit is imposed, the result is closer to AMY's; if not, 
the result is closer to XG's. This raises the concern whether this approximation 
is good for computing t]. 

This issue has been settled in Ref. [71] by using the exact amplitude for the 
23 process, which removes both the soft gluon bremsstrahlung approximation and 
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the collinear approximation to the 23process. The result of Ref. [7T] shows that 
the contribution from the 23 process lies between AMY's and XG's result but 
more close to AMY's. So the 23 process is less important than the 22 one in most 
range of coupling constant. This is consistent to the perturbative approach where 
higher order processes are only perturbation to lower order processes. 



Chapter 5 

Applications of AdS / CFT duality 



Inspired by the great success of computing the ratio rj/s of a strongly coupled super 
Yang-Mills plasma by AdS/CFT duality, many work (HQ ESQ M ESI ESQ ES] 
appear on the market about applying AdS/CFT correspondence to relativistic 
hydrodynamics with Bjorken boost invariance |196j . People try to establish a 
well-defined gravity dual to relativistic fluid dynamical by AdS / CFT duality. 



In Sec. 5.1, we investigate the shear viscosity rj of strongly coupled super 
Yang-Mills (SYM) plasma in late time of hydrodynamic evolution with Bjorken 
scaling via AdS / CFT duality. We obtain the metric g^ v in a proper time dependent 
AdSs space via holographic renormalization, whose boundary condition is given 
by energy-momentum tensor of the QGP with transverse expansion or radial flow. 
With this metric we compute rj of fluids in 0+1 and 1+1 dimension without and 
with radial flow. We find the ratio rj/s = l/(4ir) in 0+1 dimension consistent with 
the KSS bound if next-to-leading terms in proper time are included in the equation 
of motion for metric perturbations. For 1+1 dimension the result is unchanged in 
the leading order of transverse rapidity |172| . 

In Sec. 5.2, we consider a string-junction holographic model of a probe baryon 
in the finite-temperature AdS background. We investigate the screening length 
for a high spin baryon. By defining the screening length as the critical separation 
of quarks, we compute the u (spin) dependence of the baryon screening length 
numerically and find that baryons with high spin dissociate more easily. Finally, 
we discuss the Regge-like relation between the angular momentum J and the total 
energy E 2 for baryons |124j . 



5.1 Shear viscosity in late time 
5.1.1 Kubo relation 

The Kubo relation can be obtained by the statistical analysis [671 [68J. Here a 
simple way given to derive the Kubo formula given based on Ref. [50j which is 
associated with the AdS / CFT duality. 

The action with a source J a {x) and its operator O a (x) is written as 

S = S + [ J a (x)O a (x) . 

J X 

The perturbation of J a {x) gives 



50 
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(O a {x)) = - I G*(x-y)J b (x), (5.1) 



y 



with the help of the linear response theory. Here G^ b is the retarded Green's 
function defined by 

iG%(x -y) = 6(x° - y°) ([O a (x), O b (y)}} . (5.2) 

It is known that the metric g^ u as a source is coupled to T^ u as an operator 
in the general theory of relativity. For the sake of simplicity, the metric g^ v is 
considered as a homogeneous one 5ij with a perturbation hij(t) <C 1. Moreover, 
hij is assumed to be traceless ha = 0. In the local rest frame, i.e. = (1, 0, 0, 0), 



the shear viscous tensor defined in Eq.(2.28) with covariant derivatives reads 



a xy = 2rjT° xy = rjd h xy . (5.3) 

Considering a perturbation in the form of the plane wave = (h ) i je' lujt with 
constant ho yields 

(a xy ) = iujt] + 0{u 2 ) . (5.4) 



Using Eq.(5.1), one obtains 

i^xy) = — J Gxy,ijhij = — J G^ y xy h xy . (5.5) 



fluid dynamical limit, 



Substituting Eq.(|5.5|) into Eq.(5.3) and taking the low frequency limit, i.e. the 

lim- / dtdxe wt 6{t) < [T xy {x),T xy {0)} > . (5.6) 

>Ooj J 



V 

Then we derive the Kubo formula 

T) = -lim -ImG^ (u, 0) , (5.7) 

where in the low wave-number limit the Fourier transform of the retarded Green 
function gives 

G xy , xy (oo,0) = -mm j dtdxe- tk ^6(t) < [T xy (t 7 x),T xy (0,0)] > . (5.8) 

5.1.2 Hydrodynamics with Bjorken boost invariance 

In heavy ion collisions, it is assumed that all particles are created in the same 
proper time after the collisions. In the laboratory frame, one will observe in the 
center region of the collisions that the particles moving fast are produced much 
earlier than those moving slowly due to the Lorentz transformation. So one can 
assume that the space rapidity 77 = | In (f^f) is equal to the momentum rapidity 

y = I In { ^p^p j ) f° r collisions along the z-direction. After some simplification, 
the following formula is obtained 

tanhy = — = v z = - , (5.9) 
Po t 
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where r = \Jt 2 — z 2 is the transverse proper time. 

Bjorken [196J suggested that all the thermal quantities should be independent 
of the rapidity y. It is convenient to use the coordinates (r, y) instead of (t, z). Us- 
ing the coordinates (r, y) and the equation of state e = 3P, the thermal quantities 
of the ideal fluid can be solved as 



x 



Xq 



(5.10) 



e,P,m = 4/3, 



where xq is the initial value for x in the proper time To, and for x 
for x = T, m = 1/3, and for x = s, n, m = 1. 

In order to describe the evolution of the fluid with transverse expansion, besides 
(r, y) the cylindrical coordinates will also be used 



(T,y,r,6) 



(5.11) 



with r and 9 are radius and azimuthal angle in transverse plane. The velocity of 
the fluid cells can be parametrized as 



1 



dr 



f v/l-i? 2 (r,r)/' 



R(r,r) 
1,0, -^-Ao 



T 



[cosh a(r, r), 0, sinha(r, r), 0] 



(5.12) 



where the total proper time Tf is given by Tf = a/t 2 — r 2 and R(t, r) is an unknown 
function which has to be determined by the evolution equations of the fluid. For 
simplicity, the total baryon number density is assumed to vanish. The energy- 
momentum tensor reads 



T, 



( — I + |e cosh 2 a 


■Tjesinhacosha 




V 



The conservation equations are 





3 ' 






— gesinhacosha 



| + gesinh a 









\ 



(5.13) 



V W T^ = d u T^ + T^T 



where T^ a are Christoffel symbols for the metric g^J = diag(— 1, r 2 , 1, r 2 ). Solving 



(0) 



Eqs. (5.14) leads to 



(5.14) 



<9 T lne 
dAnR 



2 2 cosh 2 a 
t 2 cosh 2 a + 1 



l + - + d T R 



1R/T + 2 + 2(1 - d r R) cosh 2 a 
t 2 cosh 2 a + 1 



(5.15) 



With a given initial condition R(t , r) = £r r, where £ is set to 0.05 fm 1 given by 



Ref. [T5| I197j . the numerical results of the Eqs. (5.15) are shown in Fig. 5.1 and 
are identical to the results in Ref. [T5l 1197J. The results indicate that in this case 
the radial velocity is proportional to the distance, i.e. u r ~ £r, and is observed to 
rise sharply at the early time and fall with increasing transverse proper time. The 



evolution of the energy density shown in Fig. 5.2 is found to damp in the power 
series of r 4//3 if the transverse expansion is negligible. 
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Figure 5.2: The energy density as a function of r and r. 
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5.1.3 Holographic renormalization 

Generally, the AdS metric can be written in the form of Fefferman coordinates 



ds 2 = g MN dX M dX N = ^[g^(x,z)dx^dx u + dz 2 ] , (5.16) 

where X M = z) and 'g^x, z) is the metric tensor in the 4-dimensional space. 
In Eq.( 5.16| , there is divergences at the boundary z — > 0. To cancel this diver- 



gence, i.e. to renormalize the metric at the boundary z — >■ 0, it is straightfor- 
ward to assume that g^ ~ z 2 f(x,z) where f(x, z) is a function free of singularities 
near the boundary. The authors of |199l 1200] obtain the form of f(x, z) 

^(x,^) = ^z 2 ^)(x), (5.17) 

n=0 



and prove that 'g^ is the metric for a conformal theory. In the expansion (5.17), 
9$(%) is just the metric in the 4-dimensional flat space. The term gff(x) can be 
proved to vanish in a conformal theory. The term g$(x) is associated with the 
energy-momentum tensor of the conformal theory, 

g$(x) oc (7» , (5.18) 

where the factor s-% = 1 with N c the number of colors. Higher order terms are 
given by the Einstein equation 

Rmn - ^mnR + GgMN = , (5.19) 

with the given initial condition g^J (x) and g$(x), the cosmology constant A = 
_d(d-2) _ _g jj ere R MN and R are the curvature tensor and scalar of AdS space, 
respectively. In this case, one finds that a conformal system is located at the 
boundary of an AdS space. 

5.1.4 AdS metric with radial flow 

The authors of Ref. |191^ 1192} 11931 1194] suggest that one can use the holographic 
renormalization to regard the energy-momentum tensor of a boost invariant fireball 



as a boundary condition for g^ u in Eq.(5.17). By solving the Einstein equation 



(5.19), an AdS metric with the properties of the boost invariant fireball can be 
obtained. 

In cylindrical coordinates, 'g tiv {x^z) can be cast in the form 

g^ u (x, z)dx^dx v = -Adr 2 + Br 2 dy 2 + C (dr 2 + r 2 d6 2 ) + 2Ddrdr , 

where drdr is off-diagonal element. It is too complicated to solve the AdS met- 
ric directly in the above form. A simplification is necessary. As argued in the 
last subsection, the radial velocity is negligible compared to the transverse one. 
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Therefore, one can assume a in Eq.(5.13) is small, so T^ v can be expanded in a, 





a 



( P o 



T, 



r 2 \p 



o \ 






/ o 






\ r 2 \p J 

( y o o o \ 








-Ip o\ 





-Ip o 

\ 




) 




\p 
\ j 



+ O (a 3 ) 



(5.20) 



The off-diagonal terms proportional to a only arise in the next-to-leading order 
of the expansion. Therefore, the modification from the radial flow should be also 
in the order of 0(a). Solving the Einstein equation (5.19) with the boundary 
condition given by Eq.(5.18) and (5.20) yields 

/ 



1 



gMN 



(i-«) a 

l+a 



4 (1-a) 2 
3 a l+a 



\ 










r 2 (l + a) 






4 (1-a) 2 
3 a l+a 


l + a 









r 2 (l + a) 




o\ 







(5.21) 



with a = and eo is the same as the one in Eq.(5.10). Note here a is set to 
a scaling and higher order contribution r _1 is neglected in the late time r — > oo 
limit. 

For a = 0, the metric is 



ds 2 



(1 - z 4 /z 4 H ) 2 dr 2 1 + z 4 /z 4 H 



1 + z 4 /z 



+ 



H 



-[T 2 dy 2 + dr 2 + r 2 d6 2 ] + 



dz 2 



with the horizon of the black hole 

z H = 



-1/4 



T 



1/3 



(5.22) 



(5.23) 



The standard D3 black AdS metric can be obtained from Eq.(5.22) 

dx 2 



ds 2 = 

with the replacement 



p 



z 2 



+ 



z 4 /z 4 H z 



dp 

■Z2 ' 



(5.24) 



z —7- z 



zh -> z# 



1 + 



V2 



The Hawking temperature can also be obtained 

1 y/2 



T, 



H 



TTZ H TTZ H 



(5.25) 
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5.1.5 Evolution of rj/s 

To compute the retarded Green function, the action of gravity is necessary. Gen- 
erally, the action in the AdS space reads 



r V 

'5D 



8n 2 R 3 



J d 5 x(n 5D - 2A) , (5.26) 



where the contributions from the matter term g ilv T iW in AdS space is ignored. 
Considering a perturbation g^ v = gffl + h^ v with h z/1 = h^ z = yields 



'5D 



AT2 

8^ 
8^ 

N 2 



8tt 2 
8^2 



J <£x^g-{K^ - X -g^Tl + Ag^)h^ 
J d 5 xV=g(K$ - lg$K {0) + Agf})h»» 
j d 5 xV^g(sn, u - l - v^ ( o) - lgj$6K + Ah^)h^ 

J d 5 x^(51Z^ - \g^5TZ + 4V)^ , (5.27) 



where IZq = — 20 and A = —6 have been used. For most AdS metric, the h\ = 
is decoupled with others 



N 



<5D 



•57T 



^TT 2 



J d 5 x^(5R 12 + Ah 12 )h 12 

J d 5 x^g~ f-^d^d^ + ■ ■ ■ } . (5.28) 



In + 1-dimensional case, i.e. the radial flow vanishes, the metric gives 

ds 2 = \ ! J\~ a)2 dr 2 + (1 + a)[r 2 dy 2 + dx\ + dxl]\ + ^ , (5.29) 
z 2 y 1 + a J z 2 

which is identical to Ref. |191[ 11921 1193^ I194| . The equation of motion for <fi is 
given by 

d lx W=gg^d v <j>) = Q. (5.30) 

Assuming 

<j>(r,y,z) = <p (y,T)f(z) 

T^rri exp(-iuT cosh y + ip 3 r sinh y)<j> Q {u3, Pz)f P {z) , 

and taking the limit r — > oo while keeping a as a scaling yields 

(1 ~ Q2) ^ + /p ( l-a 2 ^ 2c ° Sh2y ~ /p(1 ~ a)( ~^ Sinhy+P3CQSh?/)2 = ° ' (5,31) 

In central rapidity y ~ and static limit p% = 0, the solution of the above equation 
is 

f p (z) = (1 - z') w , (5.32) 
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with 



z — >■ z 



1/4 



/ _^ 1 fPo 
CO — >■ — = — 

2v/2 V3 



-1/4 



(5.33) 



In this case, the action (|5.34|) becomes 
hD = 



N 2 



167T 2 
1 



doc cIjQlj 



2-0 O (-W, -P 3 )0O(W,P 3 ) 



(27T) 



X 



q 2 , r x ^/p(^) 
3 j-p{ z ) 



dz 



Z = Z H 



dp 



(f>o(-u, -p)<f) {uj,p)F(u)) 



(5.34) 



where 



N 2 1-a 



16tt 5 



AT 2 



16tt 2 3 
V2N 2 



p r- 4 / 3 lim (1 - 



;po\ 3 / 4 

za; : - I — I 



3 J 



T 



(5.35) 



The retarded Green function G R can be written in the form 

G R oc F(p) + F(-p) . 

The negative modes of the function F(p) is F(—p) = F*(p). Thus, G R oc Re F{p) 
and the imaginary part of retarded Green functions in CFT via the AdS/CFT 
duality is found to be zero. The authors of Ref. [201J assume that the imaginary 
part of G R is given by 21m F(p). In this case, the retarded Green function reads 



G? 2 , 12 (u;, 0) 



-2F(u) 



^ at2 3/4 -1 



3 3 / 4 x 4tt 



The shear viscosity can be obtained via the Kubo relation (|5.7|) 

V 



pTVni 



6 3 / 4 v^F t 



(5.36) 



(5.37) 



On the other hand, from Eq.(5.25) the entropy per transverse area and rapidity 
reads 



T2 \ 1/4 



N 2 



m 

Finally, the famous KSS bound is obtained 

T) 1 

s Att 

The metric in the 1 + 1-dimensional AdS space is written as 



(5.38) 



(5.39) 



ds 2 



1 



(1-°)' 
1 + a 



dt 2 + (1 + a)(dx3dxs + dxidxi 



8 (l-a) 2 l , , dz 2 
+ -a— Xidxidt H — 



1 + a r 



(5.40) 
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where for convenience the rectangular transverse coordinates (a+a^) is used in- 
stead of cylindrical ones (r,8). The equation of motion for <p is found to be 



the same as Eq.(5.30). Using the same assumption as in + 1-dimensional case 



<P(T,y,z) = (j)o(y,r)f(z) yields 



d 2 f p 3 + 5a df p 1 + a 



dz 2 z(l — a 2 ) dz 
whose solutions are 



2a 



1 + a 3r 



uf p = 



fp 



7T(1-W 



/- 



(1 + iu/)T(l +kj 
f p (iu)' ->■ -iu/) 



H 



OtUJ 

3r 



In this case, F(oj) in Eq.(5.34) becomes 

AT 2 I -a 2 



F{u) 



167T 2 



dfp(z) 



dz 



16tt 2 3' 



dz' 



. V^iV 2 / Po x3/4 

where the derivative of is given by 



/ -^\i3a)'/42— iw' 



;i + *u/)r(i + *u/) (1 /)2 ^3r 



aw 



-2(i-/r' +i ^ ,+2 (— )' 



-1 Z f_l)^'/4 2 - 



■2-iuo' 



[I + iu')F(l + iu') 



Finally, the ratio r]/s turns out to be 

rj 1 

S 4:71 



+ 0(a 2 ) . 



(5.41) 



(5.42) 



(5.43) 



(5.44) 



(5.45) 



5.1.6 Discussion 

We have derived a time dependent metric dual to fluid in 1+1 dimension with 
Bjorken scaling and radial flow in late time via holographic renormalization. We 
assume the transverse expansion is small and can be treated as a perturbation. In 
that case, we solve the Einstein equation and obtain the new AdS metric. With this 
metric we calculate the ratio rj/s in late time limit. We have shown that the ratio 
for fluids in 1+1 dimension is the same as in 0+1 dimension in the leading order 
of transverse rapidity. In 1+1 dimension one can introduce the shear viscosity of 
the next-to-leading order in the stress tensor in late time solution [9,11,12]. We 
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found that the correction to the shear viscosity is of the next-to-leading order as 
shown in Eq.(5.45). Compared Eq.(5.37) with rj ~ e/ m fp ~ g~ 2 T 3 ,we see that the 
coupling is as strong at the beginning as at the late time, i.e. the evolution does 
not influence the strength of the interaction. 

This technique can also be used to investigate the evolution of QGP at the very 
early time with the the scale a ~ zr + ' s ', see e.g. Ref. P02J. The local thermal 
equilibrium for the fireball has not been established because of the anisotropic 
evolution, i.e. the temperature which is a varying slowly function of the coordi- 
nates is not well-defined. It will be of interest to investigate the phenomena near 
the phase transition, e.g. evaluation of the time for the local thermal equilibrium 



5.2 High spin baryon in hot plasma 

As mentioned in Sec. |1.3.2| the Wilson loops via AdS/CFT can be used to study 
the potential between quarks and anti-quarks. Besides these, the velocity depen- 
dence of the screening length in the QGP can also be learned about in a boost 
AdS metric [203J. In Ref. |124] . the rotation dependence of the screening length 
for the baryons in the QGP is considered. The physical picture for Wilson loops 
is similar as in Ref. |124| . They differs on the end points of the open strings. In 
Wilson loops, the two end points of open strings stay in the same brane while in 
Ref. |124 j they stay in a static brane as a boundary and a probe one (also see Fig. 



5.3). The baryons live in the boundary of the AdS space. And another D5 brane 



as a probe is in the bulk of the AdS space as shown in Fig. 5.3 



5.2.1 Setup 

The AdS metric is given by 

2 2 i 

ds 2 = -f(r)dt 2 + r —dx\ + (dp 2 + p 2 d6 2 ) + -j^dr 2 + R 2 dn 2 5 , (5.46) 
R z R A jyr) 

where p and 9 are in the X\ — %2 plane. Mapping the coordinates of the strings to 
the that of the space-time 

r = t , a = r , 9 = cot , p = p{r) , (5.47) 

and using the Nambu-Goto action 

S string = ^7 J dodr a/- det [hob] , (5.48) 

with h a b = g^ ^ag^l and l/(27ra; / ) the string tension becomes 



T r VA / / r 2 \ f 1 

* / 1 / / ' 00 f / \ \ I A 



2 



s, irm , = _ ^ *y _ y-^ - f W j ^_ + j , ,5.49) 

where T is the total time and r\, r e are the cut-off and end points of the strings, 
respectively. The total action is given by 

N c 

Stotal = & string + 1 
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where the action for the probe D5 brane is 

V(r e )TV 5 



S 



D5 



(5.50) 



(27r) 5 a /3 

with V5 the volume of the compact brane and V(r e ) = y/—goo the potential for the 
brane at r = r e . 

On the other hand, a constraint condition, the so-called force balance condition 
(FBC), in the r-direction is also considered 



5> 



(0 



i=l 



with H the Hamilton of the system and 

T dr e 

5.2.2 lu dependence 



(5.51) 



(5.52) 



From Eq.(5.49), the equation of motion for p(r) is obtained 

d 



d_d_ 

dp dr dp' 



C = 



(5.53) 



where C is the Lagrangian of the action (5.49). 



condition for the numerical calculations 



Substituting Eq.(|5.49|) into Eq.(|5.51|) we can obtain the FBC as a boundary 

(5.54) 



P ( r e) 



with 




With a given value of r e , Eq.(5.53) can be solved with the boundary condition 



(5.54). The numerical result is shown in Fig. 5.4 The separation distance of 
quarks with a given r e is 



L 



''A 



p'(r)dr = 2p(r^) 



(5.55) 



The result of l q as a function of r e is shown in Fig. 5.5 The maximum value of 
l q is defined as the screening length l s of the baryon l s as shown in Fig. 5.6 We 
obtain the form of l s as 



I, 



bcu + c 



d 



where it is drawn in Fig. 5.6 with following parameters, 

0.51 



-0.85 . 



(5.56) 



(5.57) 



0.14a; + 0.22 

In hot QGP, the screening length depends on the temperature and u or spin. The 
smaller screening length of a baryon with higher spin is easier to dissociate. The 
linear velocity of quarks v = l q u>/2 is similar to Ref. |123j if a drag force in the x% 
direction orthogonal to the rotating plane is added. 
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Figure 5.3: String and a brane as the boundary represent a baryon. 
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Figure 5.4: Embedding function p{r) as a function of u with fixed r e (upper plane) 
and l q with r\ = 100 (lower plane). 
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Figure 5.5: Separation distance of quarks as a function of r e with r\ = 100. 
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Figure 5.6: Screening length as a function of spin of baryons u (upper plane) and 
linear velocity v (lower plane). 
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5.2.3 Regge behavior 

At the time when the classic string theory was considered as a candidate for the 
theory of the strong interaction, it was found that the angular momentum of 
baryons J is proportional to the square of energy E 2 as follows 

Joe E 2 , 

which is called the Regge behavior. In this model, the angular momentum and 
energy of the strings read 



J, 



string 



dC 



2na' 



dr — = , (5.5£ 



(£pV-/(r)) Uy + ^ 2 (r) 



and 



string ^~ 



dco 



2na' 



dr y J = . (5.59) 



(gpV - f(r)) (fa + ^p' 2 (r 
The total angular momentum and energy of the system are given by 

-^total N c E s t r i n g -\- Eb rane , 

Jtotal N c J s t r i n g -f- Jbrane j (5.60) 

where 

E ^e = ■ ( 5 - 61 ) 

The numerical results for J total an d E to tai as functions of u with different cutoff 



r\ are shown in the upper and lower plane of Fig. 5.7, respectively. The relation 
between Jtotai and E 2 otoal with the same u (i.e. fixed r e ) is shown in Fig. 5.8 In 
the case of the same separation distance of quarks, i.e. the fixed value of l q , the 
relation between Jtotai and E 2 otoal is shown in Fig. 5.9 and 5.10 



Finally, another physical boundary condition is considered. In order to pick 
up some configurations for baryons with the same constituents but different oj or 
spin, the configurations satisfying the following boundary condition are chosen, 



which gives 



^ = 0, (5.62) 



i'> 2 =(l-f). (5-63) 



or 

p'(r A ) = . (5.64) 
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Figure 5.8: The relation between Jtotai and E to tai with fixed r e 
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Figure 5.10: The relation between Jtotai and E tota [ with fixed l q 
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120 



Figure 5.11: Embedding function p(r) with orthogonal boundary condition (5.64) 



It is observed that the end points of strings move with the speed of light on the 



cutoff brane with the first condition (5.63), while the string is orthogonal to this 



brane with the condition (5.64). Unfortunately, it is not guaranteed to find out 



the points to fulfill the boundary condition (5.63) for small u>. The configurations 



for the strings fulfilling with condition (5.64) are shown in Fig. 5.11 



In this case, the FBC in the p-direction becomes 

dEj 



(5.65) 



with m q = ( r A —r o) the mass of quarks and Ej the interaction potential defined 

as 



N c E s t r i n g Eq -\- E\ rrane 



where 



E„ 



2-na' 



'A 



dr 



(5.66) 



(5.67) 



is the energy of free quarks. The u dependence of angular momentum and energy 



with the FBC (5.65) is shown in Fig. 5.12 The relation between Jtotai and E } 



with the FBC (5.65) is shown in Fig. 5.13 



J totoal 



5.2A Conclusion 

We calculate the u> dependence of the baryon screening length in strongly coupled 
hot plasma via AdS/CFT duality. In our model, we consider a baryon in the bulk 
space as a probe and a baryon at the boundary with high spin using the rotating 
strings. We obtain the u> dependence of the embedding function of these strings 
and their screening length. We also investigate the relation between total angular 
momentum and energy of baryons and show the numerical results in three different 



conditions. We find that these solutions with orthogonal boundary condition ( 5.65 ) 
are the best candidates for baryons with same constituents but different spins. 
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Figure 5.13: Angular momentum and energy with the FBC (5.65) as functions of 



Chapter 6 

Fluid dynamics with triangle 
anomaly 



The RHIC data for collective flows have been well described by the ideal and 
dissipative hydrodynamics gZl EH DM HH HEl [2061 EQU [208]. The correspon- 
dence of relativistic hydrodynamics to charged black-branes was investigated by 
AdS/CFT duality [146, I145j . A new term associated with the axial anomalies 
was found in the first order dissipative hydrodynamics (see, e.g., Ref. |147| about 
holographic hydrodynamics with multiple/non-Abelian symmetries, or Ref. |148| 
in Sakai-Sugimoto model). Recently the new term has been derived in hydro- 
dynamics with a triangle anomaly |149j . A similar result was also obtained in 
microscopic theory of the superfluid |15U| . This problem is closely related to the 
so-called Chiral Magnetic Effect (CME) in heavy ion collisions [ToTl IT521 [T531 ITM] . 
When two energetic nuclei pass each other a strong magnetic field up to 10 18 G is 
formed, which breaks local parity via axial anomaly. This effect may be observed 
through charge separation. Hydrodynamics in an external background field can be 
used to pin down the CME in real time simulation. However the anomalous term 
in the charge current breaks the second law of thermodynamics unless new terms 
of vorticity and magnetic field are introduced in the charge and entropy currents 

In this chapter, we try to provide a consistent description of the kinetic equation 
with a triangle anomaly. We will derive the kinetic equation to the next to leading 
order as well as the leading order correction to the particle distribution function 
arising from anomaly. These results are compatible with the entropy principle of 
the second law of thermodynamics and the charge/energy-momentum conservation 
equations. Most of the contents in this chapter are taken from Ref. [209J 

6.1 Constraining distribution function with anomaly 
compatible to second law of thermodynamics 

In this and the next sections we will consider the most simple case with one 
charge and one particle species (without anti-particles). The relativistic Boltz- 
mann equation for the on-shell phase-space distribution f(x,p) in a background 



electromagnetic field is given by Eq.(4.1). Note that C[f] contains a normal 



collision term Co[f] and a source term from anomaly Ca[/], C[f) = Co[/] + Ca[/]- 
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We assume that Ca[J] is at most of the first order, a small quantity. The necessity 
for the source term is to make the charge conservation equation hold, 

d^f = -CE^B^ = -CE ■ B. (6.1) 

Here is the charge current and E^ = u u F^ u and B^ = |e ^ ua pu v 'F al3 are electric 
and magnetic field vectors respectively, where is the fluid velocity and e^ va p = 
_ t ^uaf3 _ _^ i f or orc [ er of Lorentz indices (fivafl) is an even/odd permutation 
of (0123). However, the presence of the source term should not influence the energy 
momentum conservation, 

One can verify that the equilibrium solution of the distribution function, 

/0 = exp[/ftv(p" - QF^x v ) - 0QM -e (6 ' 3) 



satisfies the collisionless Boltzmann equation (4.1) in an external field for con- 
stant /3 = 1/T (T is the local temperature), and /io (local chemical potential 
without electromagnetic field). Here e = 0, ±1 for Boltzmann, Bose and Fermi 
distributions respectively. When (3, w M and /io are not constants but functions 



of space-time, the Boltzmann equation (4.1) is not satisfied automatically. Note 
that we can absorb —QxyUpF^ = Qx ■ E into /i so that f has the form of an 
equilibrium distribution function, 

where \i = /i — x ■ E. 

We assume that the distribution function / in presence of an anomaly is a 



solution of the Boltzmann equation with collision terms in Eq. (4.1), where f3, 
Ufj, and jj, are functions of space-time. Generally f(x,p) can be written in the 
following form, 

f( X >P) = e (u. P - Qtl )/T +X ( X , P ) _ e = fofaP) + /ifog). ( 6 - 5 ) 



where fo(x,p) is given in Eq. (6.4) and fi(x,p) is the first order deviation from it, 



fi{x,p) = -fo(x,p)[l + ef (x,p)]x(x,p). (6.6) 

It is known that a magnetic field is closely related to a charge rotation characterized 
by vorticity. So we introduce into the distribution function terms associated with 
the vorticity-induced current = \e^ va pu u d a u^ and the magnetic field 4- vector 
Bp which are assumed to be of the first order, which provide a leading order 
correction to the particle distribution function. For simplicity we will neglect 
viscous and diffusive effects throughout the paper, then the ordinary form in the 
current scheme for x{ x iP) reads, 



X(x,p) = X(p)p-uj + X B (p)p-B, 



(6.7) 
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where X(p) and Xb{p) are functions of /i, T, and u ■ p and have mass dimension 
—2 and —3 respectively. We will show that X(p) and \b{p) must depend on 
momentum otherwise they will contradict the entropy principle from the second 
law of thermodynamics. 



Using Eq. (6.5) we can decompose the charge and entropy currents and the 
stress tensor into equilibrium values and the leading order (first order) corrections 
as j» = jg + jf , S" = Sg + S? and = Tg v + T? v with 



Jo,i( x ) = q J [ t W/o,ifop)> 
S%(x) = - f [dp]p^(f ), 



S?{x) = - J [#V(/o)/i, 
Tft(x) = | [#V/o,i(x,ri, (6.8) 



where we have defined [dp] = d g ^jjafu-p) (^s * s ^ ne degeneracy factor), ip(fo) = 
f Q ln(/ ) - e(l + e/ ) ln(l + ef ) and $(f ) = \n[f /(l + ef )] = -{u-p- Qfi)/T. 
Inserting / into the above formula, we obtain the charge and entropy currents and 
the stress tensor in equilibrium, j£ = mi 11 , Sq = sw M and Tq V = (e+ P)u^u v — P 'g^ v , 
with the energy density e, the pressure P, the particle number density n and the 



entropy density s = (e + P — nfi)/T. Using Eqs. (6.6|6.7|6.8 ), we obtain 



where 



T? u = DT{u»u u + uW) + D B T{u»B u +u"B fl ) , 

= -Eteuf + frB^ + iDuf + DBB"), (6.9) 



-QJ^ = l -Q J [dp][(p ■ uf - m 2 ]/ (l + e/ )A(p), 



£b = -Q4 s = \Q f [dp][(p ■ uf - m 2 ]/o(l + ef )\ B (p), 
D = - J -f^^j[dp][(p-uf~m 2 ](p.u)Ul + eh)X(pl 

Db = -^ = ^ J[dp][(p-u) 2 -m 2 ](p.u)fo(l + ef )X B (p). (6.10) 

On the other hand, £, D and as functions of fi and T can be determined 
by the second law of thermodynamics or the entropy principle together with Eq. 
(6.1 6.2). The entropy production rate is given by 

d„(su»-^) = ^d^-u^d^ + ^j-C^E-B, (6.11) 

where C is a constant whose sign is arbitrary. We find that d^su^ 1 + S^) cannot 
be positive definite unless we make a shift to introduce a new entropy current 5^ 
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as follows, 



su" + SI - (Du^ + D B B") = su^-^ + P B B 



T 



T 



{Pu" - \if + u x T x ») - (Du^ + D B B>- 



We will use the thermodynamic relation 

d lt (Pu»)=j$d lt p-T*d li u x 

and the identities 



u»u%B x 



d»B» - 2u"E p , 



e + P 



[nB x E x + B^P), 



(6.12) 



(6.13) 



(6.14) 



to evaluate c^S 1 ^. We have used the shorthand notation = fi/T in Eq. (6.13). 
Following the same procedure as in Ref. j!49| . we obtain 



where we have defined 



ss 



L0 f 



£ s %]i-d li D + 



^P 
e + P " 



+B 11 
+E ■ uj 
+E-B 



DTn 



- d„D B + 



A 9P 



1 „ 2nD 
ess + 2D 

T e + P 

1 ess M r^ A j_ nE>B 
T k B +C T + £ + p 



P 



5S _ DsTn 
~ e + P 



-6 



(6.15) 



(6.16) 



For the constraint c^S 1 ^ > to hold, we impose that all quantities inside the 
square brackets should vanish. We finally obtain 



D 



n -V* p- c sT ^ P 



-C 



sT/i 
e + P' 



(6.17) 



Using Eqs. ( |6.16|6.17[ ), one can verify that the values of P ss and Pg S are identical 
to Ref. [149J. The difference between our values in Eq. (6.17) and those in Ref. 
[149J arises from the fact that we do not use the Landau frame while the authors 
of Ref. |149j do. By equating Eq. (6.10) and (6.17), we obtain equations for A 
and X B , 



QJ, 



21 



4 = -DT, QJ, 



21 



31 



-D B T. 



(6.18) 
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Equation (6.18) forms a complete set of constraints for A and A#. We note that 

its, we 

(6.19) 



A and Xb must depend on momentum in general. If A and Xb are constants, we 
would obtain 



Q 



J. 



21 



DT D B T -*J 31 ' 
which contradict Eq. ( |6.17 ) from the entropy principle. 
We can expand X(p) and Xb{p) m powers of u ■ p, 

Hp) = J2 x ^ u - p) 4 . a ^(p) = E x ^ u ■ py 



i=0 



i=0 



So we obtain the following expressions 

nl ~ / j A i J i+n,l, -J n \ — / J s\ ^i+n,l) 
i=0 i=0 

for n = 2,3. Here the functions J nq are integrals defined in Ref. [461 1205] . 



(6.20) 



(6.21) 



J, 



nq 



(-1) S 



1 



d 3 p 



(2g + l)!! y (2tt) 3 (m-p) 



[(« " V? ~ m 2 Y(u ■ v) n - 2q kH + ef ), (6.22) 



Using Eqs. ( |6.20|6.2l[ ) in Eq. ( |6.18[ ), we can constrain the coefficients Aj and Af . 
If we expand both A(p) and Xb{p) to the first power of u ■ p, we can completely 
fix the coefficients Ao,i and A^ from Eq. (6.18) since we have two equations for 
Aq.i and two for A^ 1; 



QJ 2 i QJ31 
J31 J 41 



A 
Ai 



whose solutions to Ao,i are 

Ao 
Ai 



<5(^21^41 — ^fl) 



-DT 



-CJ41 + DTQJ 31 
£J 3 i - DTQJ 21 



(6.23) 



(6.24) 



The equations and solutions for Aq 1 are in the same form as Eqs. ( 6.23|6.24 ) with 
replacements A ,i — > Xq 1} £ — > £b and D ^ D B - 

In massless limit, the integrals of J n k are well defined as the Poly logarithm 
functions Li n (z), 



- LL(—z) 



r s) 



fS-l 



e t /z + 1 



dt. execpt for z < — 1 



(6.25) 



When n is small, the Poly logarithm functions can be expanded in ^, so the leading 
order contribution of the solution (6.24) is then 



Aq ~ — CG\ 



rp4 



Ax « CG 2 



f5 



X B « 



A 



Af « 



Ai 



(6.26) 



where G\ and G 2 are two constants, Gi = 6075007r 2 £(5) / (d g Go) and G 2 = 12607T 6 / (d g Go) 
with G = 455625C(3)C(5) - 49tt 8 . We notice that the D terms in Eq. ( |6.24| are 
negligible, so the solutions are proportional to £. 
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6.2 One charge with particle/anti-particle 



We now consider one charge case but add anti-particles to the system. We will 
calculate X(p) and Xb{p)- Since there are particles and anti-particles, we recover 
the index Q = ±1 in particle distribution function, / — > f® and /o,i — ► fox- in 
Eq. (6.8), summations over Q should be added. In Eq. ( 6.10[ ) we also have to add 
the index Q to J^i and J*f : — » Jn\ an d Jni ~~ >" Jni ( n = 2,3), an d ad d 
summations over Q into the formula of £, £g, /}, Dg- Inserting into Eq. |6.8 ), 
we obtain the charge and entropy currents and the stress tensor in equilibrium, 

[e + P)u^u v — Pg^ v ', with the particle number 



Jo 



nit", S% 



su M and 



density n = ^qQuq, the energy density e = X]g e Q5 the pressure P = J^q-Pq, 
and the entropy density s = sq = (e+P—rifi) /T. The solutions to £, -D, .Db 
are the same as in Eq. ( 6.17[ ). Then Eq. (6.18) is modified to 



7 a,+ 

J 21 

J; 



Jo 



'21 



21 



J, 



21 



J 31 
= ~Zb, J; 



J. 

A S ,+ 



A,+ 
31 



-DT, 



31 



I 1 A B 

+ J 31 



A s - 



PbT. 



(6.27) 



Eq. (6.21) now becomes 

tA,Q _ , j\ B ,Q xB T Q 

nl ~ / j A i J i+n,li J nl ~ / , A i J i+n,li 



i=0 



(6.28) 



In the case of minimal number of coefficients we can completely fix the coefficients 
Aq,i by solving following system of equations, 



8J 2 i 8J 3 i 
aJ 3 i ffJ 4 i 



A 
Ai 



-DT 



(6.29) 



where we have used the shorthand notation, 5J n i = Jni~J n i an d vJni = Jni+Jni- 
The solutions to Aq i are 



A 
Ai 



I ( -£(aJ 41 ) + DT(5J 31 ) 
A V faJsi) ~ DT(6J 21 ) 



(6.30) 



where A = (5J2i){<jJai) — (0^31) (<^3i)- The equations and solutions for Ajj x are 
in the same form as Eqs. (6.29|6.30) with replacements Aq,i — > Ajf 1; £ — > £b and 



D^D B . 

Using integration by parts in Appendix B, 5J n k and aJ n k for n = 2,3,4 and 
k — 1 are given by 



dgT^ ( ft 3 



SJ31 
SJ 41 



Ad g T 5 

7T 2 



d g T b 



r.3 . ^ 



Li 4 (— e 

r,5 



")] 



7vr 2 

^+9^+6^ 



aJ 2 i 
cr^/31 
crJ 4 i 



7T- 



7tt 2 . /i 2 



90 + 3 + 6tt 2 



20d g T 6 



Li 5 (-e 
(6.31) 
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For massless fermions and small //, we have, 

7 4 

aJ 21 « 9C(3)G, aJ 3 i « -— — GT, aJ A1 « 225C(5)GT 2 , 

15 

7tt 4 

5J 2 i ~ 7T 2 /IG, 5J 3 i ~ 36C(3)/ZGT, 5J 4 i ~ —JlGT 2 , 

O 

7T 2 

A « —G 2 G Q T 2 JL. (6.32) 



where G = and Go = — 847r 2 £(3) + 1125^(5). The solutions have very simple 
form, 



A ~ CGij^ , Ai rj -GG 2 ^j , 

A* « ^ , Af « ^ , (6.33) 

where we have used two constants, Gi = 675 ?V' ( ' 5 ' > , Go = -r??-. We notice that 

' 1 d g Go ' z rf 9 G 



the D terms in Eq. (J6.30 ) are negligible, so the solutions are proportional to 

nfi 

_ 4 



£. Note that the quantity j^p ~ /i is also small and we have dropped it, since 
nfi « AT 4 /! 2 and e + P = « cL^T 4 . 



6.3 With two charges and particle/antiparticle 

As an example for the case of two charges, we consider adding to the system the 
chirality or an axial U(l) charge to particles. Then there are two currents, one 
for each chirality, or equivalently, for the U(1)/Ua{1) charge. For simplicity we 
assume that there is an anomaly for the axial charge current but no anomaly for the 
charge one. There are distribution functions for right-hand and left-hand particles, 
fa(x, p) (a = R,L), with chemical potentials Hr.l — H^Ha- As an extension to Eq. 



(6.7), the corrections Xa(x,p) in fa( x ,p) are now Xa(x,p) = A a p M w M + \ a BP >J 'B i 
Instead of right-hand and left-hand quantities X a , we can equivalently use X = 
Xr+Xl and Xa = Xr—Xl, where X = A, A#, £, £b> n i s i e > Pi3n- The distribution 
functions f®(x,p) satisfy two separate Boltzmann equations of the following form, 

The L r (l)/L r A(l) charge and entropy currents and the stress tensor in equilibrium 
are give by: j'q = ra/ 1 , = hau 11 , Sq = su^ and Tq V = (e + P)u^u u — Pg^", 
with the particle number density n = Y2 a Q Q n aQ (a, = R, L; Q = ±1) , the energy 
density e = J2aQ e aQ> trie pressure P = J2 aQ P a Q, and the entropy density s = 

Y,aQ S aQ = (t + P-Y / a n aVa)/T 



Similar to Eqs. (6.48 16.50 ), the Z7(1)/C/a(1) charge conservation equations (6.1 ) 
can be derived as 

d,f{x) = J[dp](C R , + -C Rt+ + C Lt+ -C Lt+ ) = 0, (6.35) 
d^x) = J[dp](C Rt+ -C R ,--C L ,+ + C L , + ) = -CE' t B lt , (6.36) 
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where we have used Eq. (6.49). The energy and momentum conservation equation 
reads, 

dpT*" - F"% = J[dp)p u (C R , + +C R> -+C L , + +C L ,-) = 0, (6.37) 



where we have used Eq. (6.52). 



Similar to Eq. (6.10), we can express £ a , £ a #, D and Db in terms of A a and 



Kb as, 



A - + -DT = J 31 



J 31 J 31 J 31 > 



for a = R,L. We have the following first order corrections, 



a=i?,L 

T^ u = DT{vPu? + u v u») + D B T(u fl B v + u p B»), 



a=R,L 



i=null,A 



- Y ^{^ + i iB B») + {D^ + D B B»). 



(6.39) 



It can be verified that the entropy current in Eq. (6.39) cannot satisfy duS^ > 



unless C — 0. In order to ensure the positivity of d^S^ in presence of an anomaly, 
one would have to substract the vector = Duj^ 1 + D B B^ from S 74 . With £7(1) 
and Ua(X) charges, we have 



a=R,L 

= ^{Pu»- £ ^ a +u x T^)-Q». 



(6.40) 



a=R,L 



In the same way as in Sec. |6.1[ the divergence of the entropy current can be 
evaluated as, 



d ^a ( 

.a=R,L ^ 



n<TD \ 2D 

£ a - duD H -<9„P 



P 



P 



+P ft 



S ^ ( ^^^P - J - 9f M D B + n "..< P 



,a=R,L 



e + P 



e + P 



-E ■ oj 



+E-B 



E 

a=R,L 

E 

_a=R,L 



n a D £ a \ 2nD 

' + — ^ - 2P B 



e + P T 7 e + P 



T e + P 



(6.41) 
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following Eqs. (6.35 (6.37)). By imposing all quantities inside the square brackets 
to vanish we can solve £ a , £ a g, D, Db as follows, 



D 



-C^ A T 
-IChua ( 1 



~ e + P 
nil 



e + P 



D B 

u = 



-Cli a li 
-Cll 2 



t 



AB 



^ _ 2ti a lia \ 
e + P J 



-Cli 1 



e + PJ 



(6.42) 



In small lia limit, Eq.(6.42) will be 



D = 
D B 



-Cll a T , 

—CflAp, 



-IClilia 
-CfU 2 , 



-Clla , 



(6.43) 



which is identical to the result of Ref. |152[ I154| . Here we have assumed that 
all integral constants are vanishing. Note that Eq. (6.43) is the result of the 



entropy principle in the hydrodynamic approach which was also obtained in Ref. 
|2imi2TTll2I2] . 

Equivalently we can use £r = (£ + £a)/2 and £l = (£ — £a)/2 to determine A a , 
X a B (a = R, L) via solving a system of equations (6.38), where the first/second line 



(each has three equations) is for Xr^l/ Xrb,lb- For minimal number of coefficients 
we can determine the values of these coefficients completely, A^ and Xrb,lb can 
be expanded to the zeroth or first power of u ■ p. For example, if we expand Xr to 
the zeroth power, then we have to expand Al to the first power, and vice versa. 
Suppose we take the former case, \r = Xro and Al = Xlo + Xlx(u ■ p), we can 
solve Xr l as 



Aro 
Alo 
Ali 



<R 



6J« 



{SJl 



1 

A 



31h 

M*J&) ~ X ( DT — i R "-f^ ) (6J&), 



4- 



'■8J&J 



(6.44) 



where A = ((^/^(crJ^) — (5 J^) (cr J^) . The solutions to \rb,lb take the same 
form as above with replacements A a — > X a s, £ a — > £,aB and D — > Db if we assume 
the same expansion as X RiL : X RB = Xrb,o and X LB = X LB fl + X LB ,i(u ■ p). 

For massless fermions and small JIrl ( or equivalently small Jt and Ji^), we 
obtain 



jL.R 

aJ 21 



SJ L ' R 



'21 



A 



7T 



: 9C(3)G , 

^ 2 p-L,Rp , 

-G 2 GqT 2 ]1l 



T L,R 



7ir 4 



GT 



15 

36C{3)p L>R GT 



a J, 



L,R 
41 



X J L > R ~ 



225C(5)GT 2 
—^-^l^rGT 2 



(6.45) 



where G and G are the same as in the former section. Then the solutions to Xr^ 
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arc 



3C fj, ix + 2/i A 3C fi ( a + 1\i A — 
Xro 1-7^ ; , X L0 — —T7^ <^4 : h Cz 5 - 



d g T 3 + fi A ' d g T 3 \ H + ^ — 

14C /i /i/iA . 3C 1 

3C 1 307r 2 C /i/M 

A M>0 ^^, A ^"^G7T^' (6 ' 46) 

where G4 = — ^847r 2 C(3) and G 5 = ^1125^(5). The coefficient ratios of X a i/X a B,i 
(i = 0,1) are proportional to /i times dimensionless factors, 

Xro A* + 2/^4 A L0 / /i + 2/i A /j,-2/j,j 
A 4 > 1 ~ A* ^4 h G 5 - 



Xrb,o A 1 + A 1 a ^lb,o V A 1 + A 1 a ^ — Ha 

ALB,1 10 A 1 — A 4 A 

Note that Alb,i <C TXlb,o ~ TXrb,o, so both Xbr{p) and Xbl(p) can be 
constants at small /Z and /Z^ limit. This property is quite different from the one- 
charge case in which Xb{p) must have momentum dependence in order to comply 
with the entropy principle. 



6.4 Collision and anomalous source terms 



In this section we will show that a general form of X(p) and Xb{p) are compatible to 
the charge and energy-momentum conservation equations ( |6.1 ) and (6.2). We will 
also derive equations for the collision and anomalous source terms. For simplicity 
we consider the single charge case without anti-particles. 

The charge conservation equation (6.1) can be derived from the Boltzmann 

d 3 p 



equation (4.1 ) as 



(2n) 3 E } 
d 3 p 



(2n) 3 E p 
where have used the identity 

d 3 p 

2 En 



-p^d^fixiP) 



d 3 p 

WFe p 



1 " u d Pv 



pfF. 



_d_ 
dpv 



f 



J d A pB{jP)5{p 2 -m 2 )p^F^f 
: - J d 4 p9(p°)5(p 2 - m^F^f 

-J d 4 pf[e(p°)5(p 2 -m 2 )F^5, u 

+ 9(p )F^2g a p P a 8%5'(p) 
+ 5(p°)F°^5(p 2 - m 2 )] 



=0 . 



d 3 p 
(27r) 3 E i 



-c[f] 



(6.48) 



(6.49) 
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Then the momentum integral of the collision term must obey 

d 3 p 



(2n) 3 E 1 



-C[f] = -CE ■ B, 



(6.50) 



so that Eq. (6.1) can hold. The energy and momentum conservation equation 

equ 

d 3 p 



(6.2) can be derived from the Boltzmann equation (4.1) as 
d 3 p 



(2tt) 3 E 1 



-p»p v d,f 

> 

d 3 p 
(2n) 3 E j 



(2n) 3 E, 



d 3 p 



dp? J (2tt) 3 E 1 



-p v C[f]., 



where we have used 



ycif] 

(6.51) 



d 3 p df 

j2^fE~ F ^ Vv dr 

d*pe(p°)6(p 2 - m 2 )F^ a PtMPl 



dp a 



= ~ J d A pf^[6(p°)5(p 2 - m 2 )F^p,p v ] 

= ~ f d 4 pf[9(p°)5(p 2 - m 2 )F» a (5, a p v + p,6 v 

+ 5(p°)F^ PlI p u 5(p 2 - m 2 ) 
+ F fia p fJ/ p u 6(p )5'(p 2 - m 2 )2g»»p p 8 va ] 
=F» a Ja . 

Then we require that the collision term must satisfy 

d 3 p 



(2n) 3 E } 



VC[f] = 0, 



so that Eq. (6.2) can hold 



We can expand C[f) to the second order as 

C[f] = C [f + f 1 + f 2 ]+C A [f]^C 1 + C 2 
where we have used the property Cq[/o] = 0, and defined 



Ci 

c 2 



dC 
df 

dCo 
df 



fi + Cai, 



f=fo 



1 d 2 C 



f=fo 



2 df 



f!+c A 2. 



f=fo 



(6.52) 



(6.53) 



(6.54) 



(6.55) 

Note that the general form for the normal part of C\ is ^ f\ — H\(u-p)p-UJ + 

f= { " 

H\ B (u-p)p- B. When inserting the distribution function (6.5) into the Boltzmann 
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equation (4.1) and using Eq. (6.54), we obtain the Boltzmann equations to the 
first and second order, 



P 



P 



d 




-F — 


dx^ 


""dp. 


d 






dx^ 





(6.56) 
(6.57) 



From Eqs. ( 6.55|6.56 ) we can determine the anomalous source term of the first 
order, 

C A i = -H x {u-p)p-u;-H XB {u-p)p-B-f {l + ef ) I fd tl [{u-p- t io)/T}. (6.58) 

By evaluating the left hand sides of Eq. (6.57), we can fix C2 as follows 

C 2 = /o(l + e/o){-[(l + 2e/ )A^'(/o) + (^A)K^ 
-[(1 + 2ef )X B d^'(f ) + (d li X B )Vp v B u 



+ 
+ 



-/?(! + 2e/ )A + 



dX 



-P(l + 2ef )X B + 



d{u ■ p) 
dX B 



p^p v E^ v + Xp^F^u" 



(6.59) 



d{u ■ p) 

-x^yd,uj u - Xbp^om ' 

Taking momentum integrals for Eqs. ( 6.56|6.57 ), we obtain the divergences of 
charge currents to the first and second order, 



d^{x) = j [dp)C x = j [dp]C A1 , 
dMx) = J[dp}C 2 . 



(6.60) 
(6.61) 



where we have used the property, 

/ [M ^ h = [ [dp] [H x {u ■p)p-co + H Xb (u ■ p)p ■ B] = 0. (6.62) 

J a J f=fo J 



With Eq. (6.8) and identities in Eq. (6.14), the left hand side of Eq. (6.61) can 



be evaluated as 



u ■ E \2£ B - 



2n£ 
e + P 
2£ 



B-E- 



e + P 
8,P ) + B» (d^B 



e + p ^' 1 - " rw;-- ^ M P). (6.63) 

In order for Eq. (6.1) to be satisfied, it is required that d^j^x) must have the 
form 



d^(x) = -CB-E-d^(x) 

2< 



-u ■ E \ 2£ B — 



e + P 
2£ 



B-E C- 



B 



e + P 



w " d ^ - "T^ V - PM dJ B - 



e + P 



^fpV)- (6-64) 
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We see that c^j'q (x) is not vanishing but of the second order though it is super- 
ficially a first order quantity. This is very important otherwise it would lead to 
9fiji{x) = —CE ■ B and give rise to incompatible results for £, £b, D, Db with 
Ref. |149| . Therefore we should keep in mind that it is the full charge current 



that satisfies charge conservation equation (6.1). We note that the expression of 



cWo( x ) m Eq. (6.64) also gives the momentum integral of the anomalous source 
term of the first order, 



[dp]C 



Al 



-u ■ E I 2£ B 



2< 



e + P 
2£ 



B-E C 



e + P 



d„P - B»l dJ B 



e + P 



P 



d^P . (6.65) 



For the energy and momentum conservation we obtain, 



[dp} P u C 2 



(6.66) 
(6.67) 



With Eqs. (6.816.14), we can evaluate the left hand side of Eq. (6.67) as 



dpT? v - F VfX jI = u u [DTd -u + u- d(DT)} + B u [D B T8 -u + u- d(D B T)} 

+«> ■ E(2D B T - ^) - E ■ + w • [d(DT) 2DT 



e + P 



P 



e + P 



+B ■ [d(D B T) - 



D B T 



e + P 



dP)} 



+DT{u -du u + u- du v ) + D B T(u ■ dB u + B ■ du v ) 

-F^^ + ^BJ, (6 

where one can verify that each term in the right hand side is of second order. The 



left hand side of Eq. (6.66) is then given by 



dX V - F U X = -(W - F^jI 



which is also a second order quantity though p u Ci 
superficially of first order. One might question the validity of Eq. (16.141) which 



rVnH l Jt F — 



(6.69) 
fo is 



follows B^Tq 1 ' — F^^j® = 0, but it is not a problem here since this equation really 
holds at the first order or the leading order but not true at the second order. It 



is essential that the energy momentum equation (6.2) hold for the full quantities 



T^ v and j^, and not for Tq\ and j®' 1 separately, otherwise the results would be 
contradictory to those of Ref. [149J following the entropy principle of the second 
law of thermodynamics. 



We now obtain the momentum integral of p u C\ t 2 from Eqs. ( 6.66|6.67 ) with 
dpTft - F^j / given by Eqs.(|6.68|6.69~|). 



6.5 Discussions and conclusions 



We have shown that induced terms related to the vorticity and magnetic field in 
the charge and entropy currents from a triangle anomaly can be derived in kinetic 
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theory by introducing correction terms to the phase space distribution function at 
the first order. We demonstrated that the anomalous source terms are necessary 
to ensure that the equations for the charge and energy-momentum conservation 
are satisfied and that the correction terms of distribution functions are compatible 
to these equations. 

As examples for the correction terms of distribution functions, we focus on 
the massless fermionic system in three cases for small fi/T, with one charge 

and one particle species (without anti-particles), with one charge 
and particles/anti-particles, and with two charges [£7(1) x C/a(1)]- In the latter 
two cases, the coefficients for u and B terms in distribution functions are found 
to be proportional to C/i/T 3 and C/T 3 respectively. In the two-charges case the 
coefficients can be constants or independent of momentum, such a property is 
impossible for the one-charge case since it is not allowed by the entropy principle. 

In the two-charges case, we assumed that there is an anomaly for the axial 
charge current but no anomaly for the charge one. The coefficients of correction 
terms for the charge/axial-charge currents and energy-momentum tensor have a 
very simple and symmetric form at small fi/T and Ha/T limit: £ w 2C/j,/j,a, 
£a « C/i 2 , £ B « Cfi A , £ab ~ Cfi, DT = -Cfi A fJ- 2 , and D B T = -C\i A li- This 
means that similar to the CME an axial anomaly can induce a residual charge 
current which is proportional to the magnetic field and the axial chemical potential 

We have a few comments about our results. In our evaluation of the correction 
terms of distribution functions, we have assumed that X(p) and Xb{p) are identical 
to particles and anti-particles. Alternatively we can assume that they have an 
opposite sign for particles to anti-particles. We can not tell which case is correct 
due to lack of deeper knowledge about these anomalous term at a microscopic 
level. Similarly we have assumed that X(p) and Xb{p) have different values for 
right-handed particles from left-handed ones. One can also assume that they 
have the same or opposite values for right-handed and left-handed particles. In 
the current framework one can not tell which is correct. Such a situation is like 
what happens in an effective theory when many effective candidates point to a 
unique microscopic theory. We also note that the solutions to X(p) and Xb{p) 
given in this paper are for the cases where the number of unknown coefficients 
in X(p) and \b(p) is equal to that of constraining equations. It is possible that 
A(p) and Xb{p) can be expanded to higher powers of (p ■ u) and then have larger 
number of unknown coefficients than that of constraining equations. In this case 
the constraining equations just provide constraints for X(p) and Xb{p) from the 
second law of thermodynamics. 



Appendix A 

Matrix elements in Eq.(|3.46j) 



For the sake of simplicity, the velocity of the fluid is parametrized as 



A X A = tanh 9 cos , 

A\\ = tanh 9 sin , 

A x 21 = (c 2 + l) sinh 2 9 cos 2 + c 2 , 

A22 = 2w sinh 9 cos , 

/la; _ AX _ 1 

^24 — ^35 — 1 ' 

^31 = ( c s + !) sin h 2 $ sin cos , 

^32 = w snin ^ sm 5 

A x 3 = w sinh 9 cos , 

= sech 2 # { sinh 4 9 cos 2 [77 + t^tt** cos(20) - t^tt** + sin(20)] 



+ sinh 2 9 [2(77 - tv^) cos 2 + 77 sin 2 0] + 77 j , 

A43 = -2^ tanh 2 # cos 2 [sinh 2 9 cos 0(71"^ cos - ir xx sin 0) + n xy ] , 
^44 = ^55 = Ttt sinh 6> cos , 



it' 



= (cosh 9, sinh # cos 0, sinh # sin 0) . 



The matrix elements of A x h are 




+w sinh 4 9 sin cos + % xy } , 
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tanh 2 #cos0 r , 9 „. __ ■> , __ , 9 , 

A 52 = , ■ , 2/3 — ^{-2smh 2 fl(7r x W0 + 27r ra sm0cos 2 

2(smh 9 cos 2 + 1) 

+n xy cos 3 0) + sinh 4 9 sin 2 (20) (tt^ cos - 2tt xx sin 0) 

-271^ Sill 0-271^008 0} , 

A x 3 = hech 2 9 {2 sinh 4 9 cos 2 0[?? - t^tt** cos(20) + T n n xx 

-T n ir xy sin(20)] + sinh 2 9[(rj + r n n xx ) cos(20) 
+3r] + T- K 'K XX - r n 7i xy sin(20)] + 2??} . 

The matrix elements of A* b are given by 

A* u = i[(c 2 + l)cosh(2^)-c 2 + l] , 

= 2 Sin ^ 9 •[ sinh 2 $ cos (2w cos 2 -|- n xx sin 2 — 7r x ^ sin cos 

(sinh 2 cos 2 0+ 1) 

+ iu sinh 4 9 cos 5 + (w + n xx ) cos + n xy sin 0} , 

7T xy COS — 7T XX Sin S 



A 13 = 2 sinh 9 ( w sin + 



A t = cos(20) 

csch 2 # + cos 2 



sinh 2 6* cos 2 0+1 



sin(20) 



csch 2 # + cos 2 ' 
^21 = (c 2 + l) sinh # cosh # cos , 
A li = (c 2 + l) sinh 6» cosh 6» sin , 
scch 3 f9 ( 

A\ 2 = — — \ 2 sinh 2 9 [(2w + tt xx ) sin 2 + 3w cos 2 - vr^ sin cos 0] 

+ w sinh 4 [cos(20) + 3] + 2w + 2tt x 
= sech 3 ${sinh 2 9 cos [iu sinh 2 9 sin + (w — 7c xx ^) sin + n xy cos 0j 

A l 2A = tanh cos , 
^25 = tanh 9 sin , 
sech 3 $ 

A\ 2 = 2 {sinh 2 ^[(w + Svr^) sin cos + 37?^ sin 2 

(sinh # cos 2 + l) 

+2ir xy cos 2 0+] + sinh 4 9[3(w + ?r xx ) sin cos 3 
+ (w + 5n xx ) sin 3 cos + 27r x?/ sin 4 + 7r XJ/ cos 4 0] 

+ sinh 6 [10 sin(20) + sin(40)l \(w - ir xx ) cos(20) 
16 

+w + 7r x:E - 7T xy sin(20)] + w sinh 8 9 sin cos 5 + 7r xy } , 
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sech 3 # 

A 33 = g 7 . r 2 fl 2 7 ~ ^ 1^ g [4(w + 20 cos(20) + (tt** - w ) cos(40) 
8 (sinh 6* cos 2 + 1) 

+21w - 97r xx + 10n xy sin(20) + tt^ sin(40)] 

+4 sinh 2 0[6w + 2^ cos(20) - An xx + 3n xy sin(20)] 

-Aw sinh 6 cos 2 [cos(20) - 3] + 8w - 8ir xx } , 

tanh 9 sin (sinh 2 sin 2 + l) 

34 sinh 2 6> cos 2 0+1 

tanh 6> cos 



2-sinh 2 fl[cos(20) -3]} , 





35 2sinh 2 #cos 2 + 2 

A\ 2 = tanh#cos0{sinh#{2sin0 [(77 - T n ir xx ) sin0 + T n ii xy cos 0] 

+77 cos 2 0} + 77 - 2r w 7r iex } , 
A43 = - tanh 6>{sinh 2 6> cos 2 [(77 - 2t 7T tv xx ) sin + 2T n n xy cos 0] 

+77 sin + 2r 7T 7r xy cos 0} , 
^44 = ^55 = ?V cosh 6> , 

A 52 = , ■ T2 a {-2sinh 2 g{sin0[-277 + r 7r 7r-cos(20) + 3r^] 

4 smh 6* cos 2 + 4 

+2^71^ cos 3 0} + sinh 4 9 sin 2 (20) [(77 - 2t^ xx ) sin + 2r^ xy cos 0] 
+4(77 - T n 7r xx ) sin - 4^^ COS 0} , 

A53 = tanh 9 j sinh 2 6> [77 cos 3 + r^7r xa: sin sin(20) - 2r 7r 7T EJ/ sin cos 2 0] 
+ (77 + T n 7r xx ) cos - T n ir xy sin j . 

The matrix elements of A y ah are 

A v n = (c 2 + l) sinh 9 cosh 9 sin , 
^21 = ( c s + l) sinh 2 sin cos , 

= . { sinh 2 # [ (w + 37r ra ) sin cos + 3tt^ sin 2 

(sinh 2 cos 2 + l) 1 

+2tt^ cos 2 0] + sinh 4 9 [3(w + tt*") sin cos 3 

+ (w + 57r xa: ) sin 3 cos + 2ir xy sin 4 + n xy cos 4 0] 

+^ sinh 6 9[10 sin(20) + sin(40)] [(w - n xx ) cos(20) 

+w + ix xx - tt^ sin(20)] + w sinh 8 # sin cos 5 + 7r X2/ j , 

A? 3 = g^^oX+l) I Sinh4 mW + 20 C ° S(20) + {nXX ~ W) C ° S(40) 
+21w - 9ir xx + 10n xy sin(20) + tt^ sin(40)] 
+4 sinh 2 #[6u> + 27r ra cos(20) - Air xx + 3tt^ sin(20)] 

-Aw sinh 6 cos 2 0[cos(20) - 3] + 8w - 8tt xx } , 
tanh 9 sin (sinh 2 sin 2 + l) 



^14 



sinh 2 9 cos 2 + 1 
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.11 tanh 9 cos r , 9 „ r . .~i 

A " = 9 • , 2fl 2 L 2 - sinh 5 [ cos ( 2 ^) - 3 } » 
2 sinh cos^ + 21- J 

^22 = ^ smn ^ sm ' 

^23 = 10 snln cos ? 
^25 = 1 5 

= (c 2 + 1) sinh 2 sin 2 + 

2 sinh # [sinh 2 9 sin cos (f)(ir xx sin — 7r XJ/ cos 0) + Tr* 1 cos + sin 

^32 = 



CI, 



(sinh 2 9 cos 2 + l)' 



Al 3 = 2 sinh 9 ( u> sin + 



7T XS/ COS ( 



sinh 2 9 cos 2 + 1 
sinh^ 9 sin 2 0+1 



^35 



sinh 2 9 cos 2 < 
sin(20) 



csch 2 # + cos 2 

= tanh 2 9 sin cos 0{sinh 2 9 [2r] + t^ix xx cos(20) - t^tt^ + r^ xy sin(20)] 
+2t7 - 2r w vr ra } , 

A\ z = — \ 2 sinh 4 9 cos 2 0[t7 + r w 7r ea: cos(20) - r^ xx + sin(20)] 

+ sinh 2 9 {??[cos(20) + 3] + 2r^ sin(20)} + 2t/ | , 

^44 = ^55 = Tvr sinh 9 SHI , 

A 52 = g / . r^-fr { sinh 2 [(t^tt- - 77) cos(40) + 9t7 + 4r.Tr- cos(20) 
8(smh 6> cos 2 + 1) I. 

-hr^ xx - 8T T 7r xy sin cos 3 0] + 2 sinh 4 9 sin 2 (20) [77 + T n ir xx cos(20) 

— T 7T 7T XX + T n 7T Xy Sin(20)] + 4 [47/ + r u 7T ra COS(20) - T n 7T XX 

-T„it xy sin(20)] + 8t7csc1i 2 # } , 
Al 3 = r n tanh 2 9 sin [sinh 2 9 sm(2(f))(iT xx sin - 71-^ cos 0) + tt^ cos - 7r XJ/ sin 0] , 

Here w = e + P. All other elements vanish. 



Appendix B 

Moments of distribution function 



B.l Basic properties 

In order to compute macroscopic quantities (energy density, pressure, etc.) via 
the microscopic distributions, the following momentum moments are of great help. 
More details can be found in Ref. [IB] with the metric g^ v = diag{ — , +, +, +} 
and Ref. [571 158] with metric same as this thesis. 
There are two kinds of moments, 



f oP a \..p a " 



' v 

J a ^- a " = f/ (lH/o)f../", (B.l) 



v 



which are called the n-th and auxiliary moments, respectively. Similar to the 
decomposition of the energy-momentum tensor T^ v with respect to an arbitrary 
velocity u^, the tensor decomposition of above moments reads, 



[n/2] 

raiQ2-«n \ r „ T jja\a 2 ...a n 

1 ~ 2-^ nq nq (<i) ' 

q=0 

[n/2] 

jaia2 ... an = Y^a nq J nq U^- a - , (B.2) 

where a nq is given by a nq = ( ™ ^ (2g — 1)!! and jj^ a ^- a n j s the rank-n projection 
operator 

jjaia 2 ...a n _ ^[a 1 a 2 ^a 2q -ia 2q u a 2q+1 yOin) 3^ 

with 



n 



n\ n(n — — k + 1) 



k J (n-k)\k\ k\ 
A{ai ...a 2qu a 2q+1 ^ u a n ) = ^\q\{n - 2q)\ 



n 

X 

permutations 



^2 A aia2 ...A a2 <*- ia2 < 1 u a2 < 1+ \..u an . (B.4) 
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These projectors possess the orthogonality property 

(2q + l)!(n - 2q)\ 



TTOt 1 a 2 ...a n TT 



TV. 



Oql ■ 



(B.5) 



Using Eq.(B.l, B.2, B.5), I nq and J nq are given by 



L nq 



J, 



nq 



(2q + l)\\J p 
■ 1 f 



\p 2 -(u- P ) 2 }«(u- P r- 2 «f A . (B.6) 



By using the identities 



nu: 



{n-2q)U^~ a ^ + 2qU {q) 

r f (»iii!...«„-2Aa„_ 1 a n _ 2 ) 
U (q) 



(?) 



(B.7) 



and Eq.(B.2), it is easy to obtain 



In+2,q — fTl Inq 

+ {2q + 3)/ n+2 

,9+1 > 

Jn+2,q = VTl J nq + (2? + 3) J n +2,g+l • (B-8) 

These integrals are Lorentz boost invariant, and therefore can be evaluated in 



local rest frame for simplicity, In that case, these integrals (B.6) will be 

And n T n+2 f 00 



1 nq 



dxx 2 ^ +l \z 2 + x 2)»/2-«-i/2( e VTO r -« _ a y 



J, 



nq 



dxx 2 ^ +l \z 2 + x 2 ) n / 2 -^ 2 - 



(2g + l)!!(27r)3y 

And g T n+2 f°° 
(2g+l)!!(27r)3y 
1 | n-2q 

o n— 1,?— 1 i ^ 1 n—\,q j 



where p 2 = m 2 with m the mass of the particles, u ■ P = E p = \/ m? + |p| and 
two new variables are defined as 



(B.9) 



m 
T 



(B.10) 



Taking variation of Eq. (B.9) gives 
dl n k 

dJ nk 



Jnkda — J n+ i^df5 , 

-g[J n -i,fe-i + (n - 2k)J n _ 1)k ]da 

-p[Jn,k-i + (n + l- 2k)J nik }d(3 , 

which are used to compute the equations of motion for the fluid. 



(B.ll) 
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The integrals of I n k and J nk can also be written in familiar form 



' nk 



And g m n+2 2 - 



kn-k 



(2g + l)!!(2vr)3 



a nM z- {k+l+1) x 



K, k+l+1 (a,z), n is even 
£ k+ i +2 (a,z), n is odd 



kn-k 



J, 



nk 



4nd g m n+2 ^ {k+l+l) 

a nklZ 



(2g + l)!!(27r)3 ^ 
£ k+ i +1 (a, z), 



n is even 



JC k+ i(a,z) + 2{k+l z +l) lC k+ i + i(a,z) } n is odd 
where the coefficients a nk \ are defined by 

(2& + 2Z + 1)!! ( \n-k 



a n ki 



{2k + 1)!! 



(B.12) 



(B.13) 



and 



K n (a,z) 



C n (a,z) 



dxx 2n (x 2 + z 2 )- 1 ' 2 



(2n- l)\\z n J e V^+? 



2n 



(2n-l)\\z n+1 J e^-"- 



(B.14) 



Similar to Eq. (B.9), /C„ is also associated with C n , 

d 



da 



d_ 
da 



C 







n+l 



-z- — (z~ n lC r . 
dz 



(B.15) 



B.2 Comparison with fluid dynamics 



We can express energy-momentum tensor and conserved charge currents of an 
ideal fluid in terms of momentum integrals 



r =/?"/ = hou" , T^= pV7 = h^u v - J 21 A 



(B.16) 



Recalling the definitions of the macroscopic quantities of fluid dynamical equa- 
tions, one obtained 

ho = n, ho = e, hx = P. (B.17) 



After using Eq.(B.9) for n = 2,3,4, more quantities can be determined 

ho + hi £ + P 



J21 - j- - j , hi 
T ho + 2/31 _ /31 



(B.18) 



From Eqs.(B.ll) the differentials of I nk and J nk are 



dn = Jioda — hod(3 
de = hoda — hod (3 



(B.19) 
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and 



n e + P 
dP = J 2 ida — Jzxdfi = —da — d/3 



which is identical to the Gibbs relation dP = sdT + fidn. 

All thermal quantities can be chosen as functions of a and 0. Alternatively, 
one can also use e and n as the thermal variables via the following identities 



da 
d(3 
dP 



D 20 
1 

1 



D 



20 



- J 20 de + J m dn) , 
'—J 10 de + -hodn) , 

[(- J21J20 + J3iJio)de + (J21J30 - J3iJ2o)dn] 



(B.20) 



where 



D n k = Jn+l,k-Jn~l,k ~ Jnk • (B.21) 

Moreover, one can also use n and s/n as thermal variables with the replacement 
of a and (3 by 



da 
dp 



D 20 
1 



D 



20 



^30 - ^^20 J ^ - J 2 onTd (- 
J21 J \n 

J20 - t^Ao ) dn - J w nTd (- 
J21 / \n 



[B.22) 



B.3 Massless limit 

In massless limit (p 2 = m 2 = 0), u -p = E p = |p|, p u p^A av = p 2 — (u- p) 2 



and Eq. (B.8) become 



1 r«; 



J, 



(2g + 3) J n>?+ i , 
(2g + 3)J n 

,9+1 ■ 



For n = 2, the equation of state for the ideal gas can be obtained 

e = 3P . 



(B.23) 



(B.24) 



Moreover, Eq. (B.9) becomes 



Ink 
J nk 



A7Td g T n+2 



n+1 



(2k + l)!!(2vr) 3 J &>-<<> - a 

dxX n+1 - ; 



A7Td g T n+2 



(2k + l)!!(2vr) 3 J ye*-* - ay 

For a Bosonic gas (a = 1), if the chemical potential is also vanishing 



(B.25) 



a = 0, 



Eq. (B.25) can be worked out 

Ink 
Jnk 



4ird T n+2 
(2^W r( " + 2K( " + 2) 



(B.26) 
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with the help of 



dxx n+l 

o e x — a 



dxx n+1 — 



: T(n + 2)((n + 2) 
- = T(n + 2)((n + l) 



(B.27) 



where T(n) is the Gamma function and £(n) is the Rieman Zeta function. 

For a Fermionic gas (a = — 1), it is necessary to consider a system with both 



particles and anti-particles. In this case, Eq. (B.25) is modified to add the con- 
tributions from anti-particles 



1 nk 



J± 

J nk 



4nd g T n+2 
(2A; + l)!!(27r)3 ./ 

And g T n+2 
(2A; + l)!!(27r)3y 



dxx 



n+l 



doc oc 



n+l 



e x ~ a - a 

e x-a 
( e x-a _ n 



i 



-1 



(B 7 28) 



where the sign (— l) n makes the differences between the charges of particles and 
anti-particles. Through the integration by parts 



n+l 



f e x-a _|_ -p 



X 



n+l 



e x-a + I 



the J^, will be related to via 



~ + rv< B+i >*" 

-a ^0 



(" + 1) pfc 

o 1 n-l,k ' 



(B.29) 



After some calculations, the integrals 7^. can be worked out 



1 nk 



And g T n+2 
{2k + 1)!!(2tt) 3 

1 



rn+2i 

m+2 1 - - 



n + l 

„ + 2 ^ V '-./ • 1 



j=0 



X V 1 " 2^ J r(2j + 2)C(2J + 2)an ~ 2 ' 
where the binomial expansions have been used 



(B.30) 
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i=o 



n 



2j + l 



The value of Jfl is obtained through Eq. (B.29). 
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B.4 Results for QGP of massless quarks and glu- 
ons 

If we neglected he mass of quarks and gluons for simplicity, in the ultra-relativistic 
QGP the integrals of moments can be worked out, 
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(B.31) 
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n2 



15 
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nO 



(B.32) 



with the help of Eq. ( |B.23[ ) where # G = iV s (JV c 2 - 1) and g Q = N s N c N f denotes 
the degrees of freedom for gluons and quarks, respectively, with N s the number of 
spin states, N c the number of color charges, and Nf the number of quark flavors. 
These results are given in Ref. [571 158] . 



B.5 Equations of motion for an ideal fluid 

As shown in Eq. ( |2.9[ ), the equation of motion of fluid dynamics are obtained from 
the conservation equations. However, the time evolution of the temperature and 
chemical potential are not known from these equations. The momentum moments 
of distribution are helpful to express the time derivative of the temperature and 
chemical potential in terms of thermodynamic quantities and fluid velocity. 
The equations of motion for an ideal fluid can be written as 

= u v d^ v 

= -J 30 p + J 31 /39 + J 20 a, (B.33) 
= A ua d u T^ 

(A v a d v p + f3u a )J 31 - {A / a d v a)J 2l , (B.34) 
= d^f 

= - J20& + J21P6 + J w a . (B.35) 
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Substituting Eq. ( |R35) into Eqs. (|R33j |R34[ ) yields 

1 







a 



D 20 
1 



D20J20 
J31 
J21 



(<^31<^20 ~~ JzqJ2iJ2o)0Q 



(B.36) 



In massless limit, i.e., J30/ J31 = J20/J21 = 3 given by Eq. (B.23), Eq. (B.36) 
becomes 



a 



-B6 , 
3^ ' 

0, 

e + P 

n 



9^(3 - + /3m a 



(B.37) 



In an external electromagnetic field F^ v = d^A u — d^A^, there is a source term 
for the energy-momentum production, 



dT p, = F uX Jx 



(B.38) 



which can also be obtained from the quantum field theory in Chap. [6] In this 
case, substituting the source term into Eqs. (B.33, B.34) yields 



- ^3o/3 + + J 20 a = 0, 

(A^/3 + /3^)J 3 i-(A^a)J 21 = nE a , 



(B.39) 



with E^ = F fJ,u u u and u ■ E — 0. The equations in present of the external field 
become 
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which in the massless limit become 
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